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On the First Boundary-value Problem 
of Linear Elastostatics 


M. E. GuRTIN & ELI STERNBERG 


1. Introduction 
The complete system of field equations in the linear equilibrium theory of 
homogeneous isotropic elastic solids, with reference to rectangular cartesian 
coordinates x; and in the usual indicial notation!, takes the form 


C15 = 3 (Mj, 5 + 4), (1) 

oO 
Tig = 2h Tnag Ci eke Cale (2) 
Tip G7 th =O § (t),= 7%). (3) 


These equations must hold throughout the region of space D occupied by the 
medium. Here u,(%), e;;(x), and 1;,;(x) are the cartesian components of the 
displacement vector-field?, the strain tensor-field, and the stress tensor-field, 
respectively, while /;(x) denotes the components of the body-force density; 
the constant parameters 4 and o designate the shear modulus and Poisson’s 
ratio, whereas 0,;; stands for the Kronecker delta. Elimination of the stresses 
and strains between (1), (2), and (3) yields the displacement equations of equilib- 
rium, which are equivalent to the single vectorial equation? 


HV ta LV EEO in D, (4) 


where V is the conventional del-operator. The present paper is concerned ex- 
clusively with the first problem of elastostatics, in which the surface displacements 
are prescribed on the boundary B of the region of space D. This problem evidently 
reduces to the determination of a vector field u(x) satisfying (4), subject to the 
boundary condition 

4 Ute) ec OULD (5) 


If D is bounded — and in the presence of appropriate regularity assumptions 
concerning B and u(x) — the uniqueness of the solution to the foregoing boundary- 


1 Throughout this paper Latin subscripts range over the integers (1, 2, 3) and 
summation over repeated subscripts is implied. Subscripts preceded by a comma 
indicate differentiation with respect to the corresponding cartesian coordinate. 

2 The single argument w represents the triplet of coordinates (%,, %2, %). 

3 Letters in boldface designate vectors. The symbols and “Xx” are used to 
indicate scalar and vector multiplication of two vectors. 
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value problem is assured by the classical uniqueness theorem due to KIRCHHOFF 
[1]4, provided the elastic constants satisfy the inequalities 


b>0,~ —ImaGa gs. (6) 
The traditional proof of KircHHorr’s theorem rests on the energy identity 


Lft-udA +i ffudv=[Wav =U, (7) 
B 


where t is the surface-traction vector with components 4;=7;;7”;, m being the 
outer unit normal of B, and U is the total strain energy. The strain-energy 
density W is given by 


1 
W = F056; = 4 (Se 26 45 Cig Sy éii), (8) 


and (6) are necessary and sufficient that W be a positive definite function of the 
components of strain. In view of (1), (7), (8), U may be regarded as a functional 
of the displacement field and admits the represention 


Ufuy=£ f (ou, Uj p+ UG 5 Mig + UG, % av. (9) 


The uniqueness theorem just cited may also be inferred from the principle 
of minimum potential energy to which it is closely related. This principle®, 
within the context of the first boundary-value problem, asserts that among 
the displacement fields satisfying the boundary condition (5) one which also 
meets the field equation (4) is characterized by an absolute minimum of the 
functional 


P{u} = Ufu} — ff-udv. (10) 


The conventional proof of the foregoing minimum principle depends once again 
on an appeal to the positive definiteness of the strain-energy density and thus 
requires (6) to hold; further, it involves the usually tacit assumption that D is 
bounded. 

It is natural to ask whether the positive definiteness of the strain-energy 
density is necessary for the truth of the uniqueness theorem or whether the 
inequalities (6) may be relaxed without loss in uniqueness. This question, though 
inconsequential from the point of view of applications of the theory to actual 
materials, is of obvious theoretical interest and, in addition, is relevant to certain 
considerations in nonlinear elasticity theory. Evidently, the requirement u>o 
in (6) may be replaced by «=-0 since KrRCHHOFF’s argument remains valid if 
W is negative-definite, rather than positive-definite. Next, it follows from 
general results due to BROWDER [3] and Morrey [4] on the uniqueness of the 
solution to the Dirichlet problem for strongly elliptic systems of partial differ- 
ential equations that (4) and (5) admit at most one (sufficiently regular) solution if 


POL le comtiont b, cteemi<ctesl, (11) 


4 See Love [2], p. 170. Numbers in brackets refer to the list of publications at 
the end of the paper. 


PRS CCE OVEN cil Dede 
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Uniqueness clearly fails when ~=0. The same is true when o=1, as observed 
by EricksEn & Tovupin [5]. Finally, it is apparent from a counter-example 
constructed by ERICKSEN [6] that }<o<1 also results in lack of uniqueness 
so that (11) cannot be further weakened ®. 

A second question arising in connection with the classical uniqueness theorem 
concerns its generalization to unbounded domains. Specifically, we confine our 
attention to regions exterior to a finite number of closed surfaces?. KIRCHHOFF’S 
argument is readily modified to accommodate such regions provided 

u;(%)=o0(1), %;(%) =O(r?) as roo, 60) be 
where r=|/x;x;. Conditions (12), however, which specify the rate of decay of 
the displacements and stresses at infinity, are highly artificial. A stronger and 
physically more useful form of the uniqueness theorem appropriate to the exterior 
first boundary-value problem was established recently by Durrin & Nott [7], 
who assume merely that u(x) tends to zero uniformly as y-+oo. Their contri- 
bution is an extension to elasticity theory of an earlier investigation by FInn & 
NOLL [8] of related issues in the theory of Stokes flows. 

The proof presented in [7] is equally applicable to bounded regions and, as 
regards limitations on the elastic constants, presupposes merely the assumptions 
(11) in place of the unnecessarily restrictive hypotheses® (6). It thus supplies 
at the same time an economical scheme for obtaining directly the generalization 
of the uniqueness theorem discussed earlier and obviates the need for invoking 
the broader, but more involved, arguments contained in [3], |4]. For this reason, 
and in order to render the present paper sensibly self-contained, it seemed desir- 
able to include in what follows, a shightly more compact version of the uniqueness 
proof due to DurFiIn & NoLt [7]. 

The proof referred to above, and spelled out in Section 2, makes use of the 
functional 


G{u} = al [x(V-u)? + (Vxu)?] dV, 
D (13) 


which plays a role analogous to that played by the strain energy U{u} in the 
conventional proof of KiRCHHOFF’s theorem. In Section 3 we deduce a minimum 
principle that is the corresponding counterpart of the classical principle of 
minimum potential energy. We show there that among the displacement fields 


6 Note that o +2 is implicit in (2) and (4); the case of the incompressible medium 
requires a separate treatment. 

7 Regions whose boundary extends to infinity are excluded from the present 
considerations. 

8 Here, as well as in the sequel, the notion of “‘order of magnitude” is used in 
its standard mathematical connotation. Thus, if v(¥) is a scalar or vector field 
defined in a neighborhood of infinity, we write v(*)=O(r*) or v(x) =0(r*) as 700 
according as |v—*v(x)| remains bounded or r—*v(%) 0 in this limit. All order-of- 
magnitude statements will henceforth be understood to refer to the limit as 700. 

9 This fact goes without comment in [7] which is chiefly concerned with the 
removal of the order-of-magnitude restrictions (12). 
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conforming to the boundary condition (5) one which in addition meets the field 
equation (4) yields an absolute minimum of 


Wiu} = G{u} — ff-udv. (14) 


This minimum principle is valid for the extended range 
pio= 0-0! = comtrord 95 in ho Kp, (15) 


of the elastic constants and applies also to the exterior problem if u(x) =0(1) 
and provided the body forces are suitably restricted’®. The Duffin-Noll theorem 
of Section 2 may be considered as a corollary of the theorem established in 
Section 3 if u>0. 

The connection between the two alternative minimum principles is explored 
in Section 4, where we show that for a bounded D the difference between the 
functionals Ufu} and G{u} depends solely on the boundary values u*(x), and 
hence is an invariant of the class of admissible displacement fields. It follows 
in particular that the principle of minimum potential energy continues to hold 
if (6) is replaced with (15), despite the fact that the strain-energy density need 
then no longer be positive definite. 

Section § contains a theorem which yields lower bounds for the functional 
W{u}. Since the structure of G{u} is simpler than that of U{u}, the results 
obtained in Sections 3, 4,5 are apt to be of interest also in connection with 
the application of direct variational methods to the first boundary-value problem 
of elastostatics. 


2. The uniqueness theorem of Duffin and Noll 


By a regular region of space D+ B we shall mean a region of space D whose 
boundary B consists of a finite number of non-intersecting “‘closed regular 
surfaces’’, the latter term being used in the sense of KELLocc [9]". Thus B 
is necessarily piecewise smooth but may have corners and edges. Further, D 
need not be simply connected or bounded. If D is infinite, however, B still 
remains bounded. The uniqueness theorem dealt with in [7] may be stated in 
the subsequent somewhat more inclusive form. 


Theorem 1. Let D+ B be a regular region of space. There exists at most one 
vector field u(x), defined and twice continuously differentiable in D+ B, which 
satisfies (4) and (5), provided (11) hold and if u(x)=0(1) in the event that D is 
infinite. 

The proof of the corresponding theorem in [7] may be adapted and condensed 
as follows. It is evidently sufficient to show that 


V2u+ aoa LY - thes Oly inp Dy (16) 


u=0 onB, (17) 


ay The generalization of the principle of minimum potential energy to exterior 
domains is elementary only if the artificial regularity conditions (12) hold and provided 
the body forces are properly behaved at infinity. 

1 See [9], p. 112. Note that the definition of a “regular region of space’’ used 
in {9] is narrower than that employed here. 
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imply u(x)=0 in D, if u(x) has the regularity properties stipulated above. 
Suppose first that D is bounded. From the divergence theorem, with the aid 
of (16) and elementary identities in vector analysis, follows the identity 


bof euV-utux(7xw)]-ndd = Giu}, (18) 


B 


where G{u} and x are given by (13). In view of (17) and (18) one has 
G{ul=0. (19) 


Since x>0 according to (11) and the second of (13), we conclude from (19) and 
the smoothness of u(x) that V- u=Vxu=0 in D, whence 


V2u=o0 inD. (20) 


But (20) and (47), in conjunction with the uniqueness theorem for the Dirichlet 
problem, assure that u(x) =0 in D. 

Next, consider the case of an infinite D and let S(o) be the surface of a sphere 
of radius g, centered at the origin and containing B wholly in its interior. On 
applying the identity (18) to the region R(e) bounded by B and S(g) and using 
(17), (13), we arrive at , 

f [xuV-utuxWVxu)|-ndA=f |x(V-u)?+ Vxu)?] dv. (21) 
S(@) R(Q) 


Operating on (16) with the divergence and the curl, respectively, one confirms 
ae V2V-u=0, V2Vxu=0 inD. (22) 
We now invoke Lemma 1 of [7], which asserts that (22) and u(x)=o0(1) imply 

V- a= bi hee (1?) (2338 


It is clear from (23), since u(x)=o(1), that the surface integral in (21) tends 
to zero as 9->oco. Consequently the same is true of the volume integral on the 
right-hand side of (21), from which the conclusion w(x) =0 in D follows as before 
in view of u(x)=o0(1). This completes the proof of Theorem 1. 


3. An associated minimum principle 


Definition 1. Let D+B be a regular region of space. Let u(x), defined and 
twice continuously differentiable in D+ B, meet the boundary condition 


u=u*(x) on B, (24) 


and let u(x)=0(1) if D is infinite. The totality of all such u(x) constitutes the 
class [tt] of kinematically admissible vector fields corresponding to D+-B and to 
prescribed surface values u*(x). 

We now state a minimum principle which is closely related to the uniqueness 
theorem of Section 2. 


12 As is well known, any solution of (16) which is twice continuously differentiable 
has continuous derivatives of all orders. In this connection see also [10]. 

13 The proof of the second of (23) is carried out in detail by Finn & Nott [8]; 
the first of (23) may be established by strictly analogous means. 
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Theorem 2. Let D+ B be a regular region of space. Let u(x) be a vector freld, 
defined and twice continuously differentiable in D+-B, which satisfies (4) and (5). 
If D is infinite, let w(x)=0(1), f(x) be continuously differentiable in D, and 


V-f=0, Vxf=0 mD. (25) 


Let [tt] be the class of kinematically admissible vector fields corresponding to D+ B 
and u*(x). Then, assuming (15) to hold, 


Tt) = 7 Pia, (26) 
where V is the functional defined in (14), and this absolute minimum is assumed 
by Y{a\ only for u(x) =u(x) im D. 

To establish this theorem it suffices to show that 


Va > Pfu}, (27) 
Pia} = {ul if and only if w(x)=u(x) in D. (28) 


Defining u’(x) through 
u'(x) = u(x) — u(x), (29) 


we consider first the case in which D is bounded. From (13), (14), and (29) 
follows directly 


Pa} — P{u} = Gf} — ff-w dV + L{u,w}, G0) 
D 
where 


Liu, u}=uf [xV-u)V-w) + (Vxu)-Vxw’)] dv. (31) 


By hypothesis and (29) both w(x) and w’(x) are twice continuously differenti- 
able in D and hence obey the elementary identities 


(V-u) (V-w') =V-(wV-u)—w-V(V-u), 
(Vxu)-Vxuw') =V- [ux Vxu)|+u-VxVxu), (32) 
VxVxu)=VV-u—V? wu. 
With the aid of (32) we may write (31) in the form 
OT Rah OF) i: V-[xu'V-u+w x(Vxu)] dv 
be (33) 
= | wu’. la u + Pare VV - ul dV. 
D 
Applying the divergence theorem to the first integral on the right-hand side 
of (33) and bearing in mind that u(x) is a solution of (4), we reach 
Liu, w} =u f [ew V-utwx (Vxu)]-ndA+ff-w dv. (34) 
B D 
Substitution from (34) into (30) yields 
Pa} — Pu} = Glu} +f ew V-ut+w x(Vxu)]-ndA. (35) 
B 
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But, by hypothesis, w(x) and &(x) assume the same boundary values u*(x), 
so that from (29) 
wu’ =0 on B. (36) 
Therefore (35) reduces to 
Pia} — P{u} = GC {u'}. (37) 


By virtue of (13) and (15), G{w’} is non-negative, whence (27) is an immediate 
consequence of (37). Further, according to (13) and (29), CU }=0 when 
a(x) =u(x). Next, G{u’}=0, in view of (13) and (15), implies V- wu’ =Vxu'=0 
in D, whence by the last of (32) 

AE A Ui hoe Oy (38) 


Invoking again the uniqueness theorem for the Dirichlet problem, we draw from 
(38), (36) that u’(x)=0 in D. Consequently, a(x) =u(x) if and only if G{u’}=0 
and thus (37) implies also (28). This completes the proof when D is bounded. 

If D is infinite, on the other hand, the identity (35) may be applied to the 
bounded region R(e) whose boundary is B+ S(@), where S(o) is once more the 
surface of a sphere of radius @, centered at the origin and containing B in its 
interior. Thus and by (36), which continues to hold in the present circumstances, 


RA amiey a Gade Lita Pak -utu xVxu)|.ndA. (39) 

Recalling (13) and (14) we note that i functionals Y and G are now given by 
{ul = G{u} — aie udV, (40) 

G{u} = 2 fle (V-u)?+ Vxu)?] dV. (41) 


In view of the assumed Cars of f(x) when D is infinite, we are entitled 
to operate on (4) with the divergence and curl, respectively, provided 7 is suf- 
ficiently large!*. Making use of (25) we find in this manner that (22) and (23) 
hold again. Also, since V-w and V xu are harmonic in a neighborhood of 
infinity and O(r?), VV-u and Vx(Vxwu) are both O(7%), so that the last of 
(32) implies that V?u=O(r-3). It now follows from (4) that 


PaO 7s), (42) 
Finally, by hypothesis and (29), 


Wo, © t=). US OT): (43) 


According to (23) and (43), the surface integral in (39) tends to zero, while 
(23), (42), (43) and (40), (41) assure that W{u} approaches a finite limit as 
o—>ce. Further, as is apparent from (40) to (43) and x>0, either G {a}, and hence 
Pia}, tends to a finite limit or W{t@}—>-+ 00 as goo. In the latter instance 
(27) is trivially met. We may therefore confine our attention to admissible fields 


we See Footnote No. 12. To demonstrate the availability of continuous third 
derivatives of u(x) in the present circumstances one merely needs to exhibit a parti- 
cular solution of (4) which has this degree of smoothness in a neighborhood of infinity. 
Such a particular solution is readily constructed with the aid of spherical harmonics 
when f(x) meets the hypotheses of Theorem 2. 
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which render G{u} in (41) convergent. For such a choice of #(x), (39) together 
with the preceding observations permit the conclusion that G {u’} is convergent 
as 9—>co, whence (37) remains valid if D is infinite. The conclusions (27), (28) 
may now be reached in precisely the same manner as before. This completes 
the proof of Theorem 2. 

Theorem 1 becomes a corollary of Theorem 2 if!® ~>0. To demonstrate this 
one need merely show that f(x)=u*(x)=0 in Theorem 2 implies u(x)=0 
throughout D. In these special circumstances u(x) =0 is an admissible vector 
field and W{u}—G{u} according to (14). Theorem 2 now implies G {u}SG {0}, 
while (13) yields G{0}=0. Hence G{u}<0. On the other hand, G{u}=0 as a 
consequence of (13). Therefore G{u}—0, which is possible only if u(x)=0 in 
D since at present u=u*(x)=0 on B. 


4. Relation to the principle of minimum potential energy 

The minimum principle asserted in Theorem 2 passes over into the classical 
principle of minimum potential energy (for the first boundary-value problem) 
if the functional Y in (26) is replaced with the potential energy ®, defined by 
(10), provided (6) are required to hold in place of (15) and on the assumption 
that D is bounded. In view of (10) and (14), the transition from ¥Y to @ is, in 
turn, effected by substituting U for G in (14). Both minimum principles are 
equivalent to the same boundary-value problem (4), (5), Equation (4) in either 
case being the Euler equation of the corresponding variational principle. This 
observation suggests a connection between the functionals G and U which we 
establish presently. 


Theorem 3. Let D+ B be a bounded regular region of space. Let u(x) and 
a(x) be defined as in Theorem 2, and set 
Hiuy = G (uy — Ue}, (44) 
where G and U are the functionals given by (13) and (9). Then 
Tava (constant), (45) 


ie. H{tt} 1s an invariant of the class of kinematically admissible vector fields [we]. 


With a view toward a proof of Theorem 3 we observe first on the basis of 
(13) that G{u} may be written as 


G{aj= 4 | (8; i, +, 0; ; — Gi, i, ) dV. (46) 
D 
From (9), (44), (46) and elementary manipulation follows 


) 
Using the divergence theorem in conjunction with (47), we are led to 
Hti}=ufhda, (48) 
B 
where ; PRs it. 
h(x) = (u; ;U; arc U; ;U;) nN; . (49) 


1° When ww <0, Theorem 2 gives w 


ay to the analogous maximum principle, from 
which Theorem 1 may be deduced. 2 : nt 
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We now complete the argument by showing that the values of h(x) on B 
depend solely on the boundary values w*(x), which are common to all members 
of [uw]. It is clearly sufficient to confirm this fact for an arbitrary regular point 
of B. Let Q be such a point, whence B has a uniquely defined tangent plane 
at Q. For convenience, choose the rectangular coordinate frame in such a way 
that the origin is at Q and the plane x,=0 coincides with the tangent plane 
of B at Q, the x,-axis pointing in the direction of the inner normal of B at Q. 
For this choice of coordinates (49) yields 


~ Ou ~ OU . (OU Ou 
Diy Vesey ee 1 Us ( bee a 
WN % ewe Te Ox, ae 0%» its Ca oak (50) 


In a neighborhood of Q the boundary B admits the parametrization 
%=€,, %, = bo, X= 9lE,,62), (51) 


where (é,, €,) is continuously differentiable and 


tod ay iz ab 
C10, 0) = Os FERN (0) ae ee) (52) 
Now ] 
U(X, %q_, Xz) = U*(G,§) on B, (53) 


in the neighborhood of Q under consideration. Also, by (51), (52), (53) 


u(x) CUS) = 5 
O72 =| (0.0)(0) pC Gann (080) Cet (54) 
Hence, from (50) and (54), 
_[,,* 04s | ,* Oug {* ( Our ouz | 
TOO) eae age = Melcaes + Tapa af¢ (55) 


Therefore h(x) at Q depends exclusively on the boundary values u* of & and is 
independent of its interior values. Recalling (48) we conclude that (45) is bound 
to hold. 

Theorem 2 is evidently an immediate consequence of Theorem 3 and of the 
principle of minimum potential energy, if D is bounded and the elastic constants 
conform to (6). Conversely, when D is bounded and with limitation to the first 
boundary-value problem, the principle of minimum potential energy follows from 
Theorem 2 and Theorem 3 even if ~ and o satisfy merely the weaker restrictions 
(15). The classical minimum principle is thus seen to be valid also when 
—oo<o<—1 or 1<0<0, in which case the strain-energy density is no longer 
positive definite. 

5. A lower-bound theorem 


The minimum principle of Theorem 2 supplies upper bounds for the functional 
Y{u} and hence, if D is bounded, also for the potential energy 


P{u} = P{u} — k, (56) 


once the constant k has been computed corresponding to a particular set of 
surface data u*(x). We now seek to determine lower bounds for Y{u}. This 
may be accomplished with the aid of a familiar scheme going back to an idea 
of Trerrtz[//] that has given rise to a large number of more general 
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investigations concerned with the estimation of quadratic functionals?®. In 
preparation for the theorem to follow we introduce the subsequent auxiliary 
definition, which is a counterpart of Definition 1. 

Definition 2. Let D+B be a regular region of space. Let u(x), defined and 
twice continuously differentiable in D+ B, meet (4), as well as the integral condition 


f [x(t — @) Va + (u* — @) x(Vx @)] nd A =0, (57) 


where x is given by the second of (13). If D ts infinite, let u(x) =o (1). The totality 
of all such iu (x) constitutes the class [tt] of statically admissible vector fields correspond- 
ing to D+ B, f(x), and to prescribed surface values u*(x). 

Theorem 4. Let D+ B, u(x), and f (x) have the same properties as in Theorem 2. 
Let [wu] be the class of statically admissible vector fields corresponding to D+ B, 
[ (x), and u*(x). Then, assuming (15) to hold, 

W{u} = max Via}, (58) 
[tt] 


where V is the functional defined in (14). 
In proving Theorem 4, one may follow closely the procedure adopted earlier 
to establish Theorem 2. At present we need to show that 


Pia} = Pu} if w(*)=u(x) in D. (60) 17 
Define w(x) through uh ne 64) 


and assume initially that D is bounded. The same reasoning that led to (35) 
now yields 
Yiu} — iw} =Glw’ 34+ uf lxeu’V-t+u"xVxi)]-ndA. (62) 


B 
Since the surface integral in (62) vanishes by virtue of (61), (57), and (5), we have 


iu} — iat = Glut, (63) 
from which (59), (60) follow at once because of (13) and (15). 
In the event that D is infinite, (63) may evidently be replaced with 


Wuj —Va@=Gw +f leu’ V-a+u’xVxa)]-ndA, (64) 
S (9) 
in which Y and G are currently defined by (40), (41), while S(o), as well as 


R(@), retain their previous meaning. 
By hypothesis and from (61), 
u=o(1), t@=o(4), w’=o(1). (65) 
Also, arguing as in the analogous part of the proof of Theorem 2, we confirm 
once again the validity of the estimates (23), (42), and in addition arrive at 
V-u=O(r*), Vxu=O(r?). (66) 


ms We shall not attempt to list here even the most significant contributions to 
the literature on this subject. Many of the pertinent references may be found at 
the end of a monograph by SyncEr [12]. 

“ Note that, in contrast to (28), the equality in (59) may hold even if a (7) = u(x) 
fails to hold throughout D. 
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According to (66) and the last of (65), the surface integral in (64) tends to zero 
as g—>oo. Similarly, using (41), (23), (61) and (66), we conclude that G{u’’} 
approaches a finite limit as @—>oo. Thus (63), and consequently (59), (60), 
continue to hold when D is infinite. The proof of Theorem 4 is now complete. 


Note added in Proof. The authors are indebted to Dr. R. J. Knops, who recently 
drew their attention to a paper by Kevin [Sir Wirt1amM TuHomson, On the reflexion 
and refraction of light, Philosophical Magazine, 26 (1888), 414]. This investigation 
is concerned in part with the relaxation of the usual conditions (6) for the stability 
of the unstrained equilibrium configuration of an isotropic elastic solid. Specifically, 
KELVIN shows that the strain energy U{u}, given in (9), coincides with the functional 
G{u}, defined in (13), provided the region occupied by the solid is finite and if the 
surface displacements vanish. He observes further that G{u} is positive provided 
the velocities of irrotational and equivoluminal waves appropriate to such a solid 
are both real, 7.e. if the inequalities (15) hold true. The preceding observations of 
KELVIN are readily extended to yield a proof of the generalized uniqueness Theorem 1 
for the case of a bounded domain. Also, the introduction of the functional G{u}, 
which in the present paper is attributed to Durrin & Nott [7], should be credited 
to KELVIN. 
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Linear Functional 
Electromagnetic Constitutive Relations and Plane Waves 
in a Hemibedral Isotropic Material 


R.A. TOUPIN & R.S. RIVLIN 


1. Introduction 
In any Lorentz frame of reference (x', t), where x* ({=1, 2, 3) are rectangular 
Cartesian spatial coordinates and ¢ is the time, MAXWELL’s equations in a 
stationary, polarizable and magnetizable medium can be written in the form 


auif a 
& [ Beda QE dl=o, (1.1) 
Bo LL 
¢ B-da=0, (1.2) 
A 
OH dl— [ D-da— [ Jp-da, (1.3) 
gL AD dD 
6D-da=f Q¢dv, (1.4) 
D=s,E+P, H=,'B—M, (485) 


where B is the magnetic flux density, E the electric field, H the magnetic intensity, 
D the electric displacement, Jp the density of free current, Qp the density of free 
charge, P the polarization, M the magnetization, and ¢) and uy are fundamental 
constants whose product & (4j=c * where c is the velocity of ight in vacuum. 
dl denotes an element of an arbitrary closed contour ¥ in space and da an 
element of a surface / bounded by ¥. 

The system of equations (1.1) to (4.5) is underdetermined, and to obtain a 
determinate system of equations amongst the twenty-two unknown field com- 
ponents (B;, £;,H;,D;, PB, M;, Jz, Qz) it is necessary to supplement the basic 
Maxwell equations with certain constitutive relations. The form of these con- 
stitutive relations depends on the nature of the material medium in which the 
electromagnetic field (H, B) resides. The great mathematical variety of these 
electromagnetic constitutive relations makes possible a great variety of physical 
phenomena embraced by and consistent with MAXWELL’s framework of equations 
(4:4) 10L(405)" 

The simplest of all media, vacuum, is characterized by the constitutive 


relations 
Jr = 0; =M=P=0. (1.6) 
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The next simplest medium is the rigid, linear, stationary, non-conducting 
dielectric, for which the constitutive relations 


py=dp=0, 
D=e-E, B=vp-H, 


were given by MAXwELL [1873, § 784]. Here € and w are constant second order 
tensors, proportional to the unit tensor if the medium be isotropic. As is known, 
constitutive relations of this simple type do not account for the observed absorp- 
tion and dispersion of electromagnetic waves in non-conductors, nor does 
MAXWELL’s device of including a linear law of conduction, J;=C - E, replacing 
(1.7)2, suffice to account for the observed magnitude of the ‘‘dielectric losses”’ 
when a material is placed in a variable electric field (Gr WHITEHEAD [1927, 
Lecture 1]). 

HopKINSON [1877], following a suggestion of MAXWELL, proposed a con- 
stitutive relation for the electric displacement in a dielectric having the form 


(Asa) 


D(i) =cE(t) +f p(¢— 2) E(x) az, (1.8) 


where the function p(w), w=0, is a decreasing function of wu. By suitable ad- 
justment of the memory function ~, HOPKINSON was able to correlate his data 
on the residual charge of Leyden jars. In proposing (1.8), HOPKINSON was 
guided by the earlier mechanical constitutive relation of BOLTZMANN [1874] in 
which the torque and twist of a wire are related by a formula identical with (1.8). 

MAxweELt (cf. RAYLEIGH [1899]; WuiTTaKER [1951]) initiated also the 
method which enjoys current favor of deriving electromagnetic constitutive 
relations from a molecular model of a material medium and the laws of mechanics. 
Elaborate theories of electromagnetic absorption and dispersion of this general 
type were developed by DRuUDE [1893, 1902], VorcT [1899], and many others. 
A detailed and modern theory of electromagnetic constitutive relations as derived 
from the dynamics of an ionic crystal lattice, including an appropriate appli- 
cation of quantum mechanics, is given in the book by Born & HUANG [1954]. 

VOLTERRA [1912] extended the Hopxkinson’s relation (1.8) to the nonlinear, 
anisotropic, and magnetic case. VOLTERRA’S most general expression of the idea 
took the form of the equations: 


D(a, t) =e€- E(z, t) + FLE(a, Tals 
“29 (1.9) 
B(x, t)=p-H (x,t) + (H(z, t)], 


where (1.9), reduces to (1.8) if the functional F is linear, isotropic, and satisfies 
certain other conditions of which we shall say more later. VOLTERRA’S theory 
was developed further by GRAFFI [1927, 1928]. 

VOLTERRA’S paper, though it treats the general anisotropic case and expresses 
a very general view, makes the traditional a priori separation of electric and 
magnetic effects characteristic of MAXWELL’s simple relations (1.7). That this 
separation is inadequate to account for numerous known phenomena such as 
optical activity and the rotation of the plane of polarization of light by a strong 
magnetic field (Faraday effect) should have been apparent from the earlier 


190 R. A. Toupin & R. S. RIVLIN: 


theory of VorcT cited above. In this paper we wish to explore some of the more 
obvious consequences of electromagnetic constitutive relations having the general 


form t 
D (x,t) =®[E (x,t), B(x, 7)], 
a (1.10) 
H (x,t) = [E (x, t), B(x, 7) ]. 


We shall show that if VoLTERRA’S relations (4.9) and the relations (1.10) are 
both linearized, then, for holohedral isotropic materials, they coincide so that 
there is, for these materials, a separation of effects; however, we shall show that 
the phenomenon of optical activity, which does not occur in holohedral isotropic 
materials, depends in an essential way upon the appearance of B in (1.10), or 
of E in (1.10). We shall show that, so far as the dispersion and absorption of 
plane electromagnetic waves is concerned, linear functional constitutive relations 
of the form (1.10) which are consistent with VOLTERRA’s principe du cycle fermé 
(explained in § 2) lead to dispersion formulae similar to those obtained by BoRN 
& Huanc from a molecular model. The simple and direct mathematical reasoning 
we apply to the relations (41.10) should have application in any theory of the 
electromagnetic field in stationary matter which leads ultimately to relations 
consistent with (1.10). The physical idea expressed by these relations is simply 
that, in a stationary medium, the values of the electric field E(x, t) and the 
magnetic flux density B(a#, t) for all times tS<# preceding the instant ¢ uniquely 
determine the values of the electric displacement D(a, ?) and the magnetic 
intensity HM (a, ¢) at the time ¢. 


2. Restrictions imposed on the constitutive relations 
by material symmetry 


If the material considered has some symmetry, the constitutive relations 
(1.10) must be form invariant under the group of transformations {S} describing 
the symmetry. For solids, the group of material symmetry transformations {S} 
is either the full orthogonal group or a subgroup of it. Let S=||s;,|| be a generic 
transformation of this group. Then the functionals on the right-hand side of 
(1.10) must satisfy the equations 

t t 


®,[E(z), B(r)] =s;; ®,[E (x), B(x)], 


—oo —oo 


tis ct (2.1) 
¥{[E(t), B(t)] = det Ss, % (E(x), B(z)], 
where ( 
E;=3;;E;, B, = det Ss,,B,, (2.2) 


for every element S of {S!. The factor det S= +44 appears in (2.1), and (2.2), 
because, as one can recall from elementary electromagnetic theory, the invariance 
of MAXWELL’s equations requires that D and E transform as absolute vectors 
under time-independent orthogonal transformations if the current and charge 


are absolute quantities, while B and H must transform as axial vectors under 
this group. 


Hemihedral Isotropic Materials in Electromagnetism 191 


In this paper, we shall treat only the case where ® and W are linear functionals. 
This will exclude the more interesting non-linear phenomena such as the Faraday 
effect, but we choose to dispense with the simpler linear theory first. In the 
literature on dielectrics (cf., ¢.g., WHITEHEAD [1927]), linearity of the functional 
® is justified by what is called the principle of superposition. Actually, we 
shall make the more restrictive assumption that ® and W are linear functionals 
having the form 


p p t t 
Dit) a eras (2) ates () + fei (t,t) E(t) de + J pelt, t)- B(x) dt, 


y= 


P p t t 
H(t) = 2.d,-E® (i) + 2b. B”() + J vale t)-B(t) dt + J po(t, t)- E(t) dt, (2.3) 


v=0 v=0 


where a,, b,, ¢,, and d, are constant tensors and the kernels @,, pz, ,, and Wy 
are continuous tensor functions of ¢ and t such that each satisfies an order 


condition of the form C 
(OPA 3 geen ee SEES ée>0. (2.4) 


(t—t)1re? 

Thus, in VOLTERRA’S terminology [1930, Chap. I], ® and W are linear functionals 
with order of continuity # having ¢ as an exceptional point. Now VOLTERRA 
[1912] has proven that, with the order condition (2.4), D(t) and H(t) will be 
periodic functions of ¢ whenever E(t) and B(t) are periodic functions (principe 
du cycle fermé) 7f and only if the kernel functions @,, ..., YY. depend on the two 
variables ¢ and t only through their difference ¢—vt. Again, we shall delimit 
the class of materials under consideration by assuming (2.4) and that D and 
H are periodic when E and B are periodic. This latter assumption rules out 
something like an “‘ageing’’ or deterioration of the material. Heredity of the 
nature we assume here is also called invariable heredity. 

For the linear functionals (2.3), the restrictions of symmetry (2.1) are satis- 
fied if and only if the second-order tensors which completely determine the 
functional relations (2.3) satisfy the restrictions: 


Sa,S7?=a,, Sb,S4=b6,, Seo,S1=q,, Sw,St=y, 
(det S) Se,S1=c,, (det S) Sd, S42 =d,, (2.5) 
(det S)Sq.S4=@s, (det S) Sp,S4=th,, 


for each element S of the group of orthogonal transformations {S}. In other 
words, a,, b,, p,, and :p, must be invariant absolute second-order tensors under 
transformations of {S} and ¢,, d,, p,, and , must be invariant axial second- 
order tensors under transformations of {S}. One can, in principle, determine 
the most general solution of the invariant theoretic problem posed by (2.5) for 
any group of orthogonal transformations by methods which are known. 

If {S} is the full orthogonal group, the material is called holohedral isotropic, 
there are no invariant axial second-order tensors, and the absolute invariant 
second-order tensors are all proportional to the unit matrix. Thus, for holohedral 
isotropic materials, the constitutive relations (2.3) must reduce to VOLTERRA’S 
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Pee p : ’ 
[1912] form D(t) = da, E(t) +S ai(t— 2) E(2) dt, 
y=0 nee 


(2.6) 
P 


H(t) = >) 8, B®) + J w(¢— 1) BO) dt. 


Thus we see that linearity and holohedral isotropy are sufficient conditions to 
eliminate B from (1.10), and E from (1.10),, as had been asswmed possible by 
VoLTERRA for a material of general symmetry. 

If {S} is the proper orthogonal group, consisting of all orthogonal trans- 
formations with positive determinant, then all of the tensors in (2.5) may be 
arbitrary multiples of the unit matrix and the constitutive relations (2.3) reduce 
to the form 


D(t)= > a, BE” (t) + : c, BY (t) + f [p(t — 1) E(t) + @2(¢ — t) B(t)] dt, 
v=0 v=0 aoe (2.7) 
H(t) = > d, B® (t) + > b, B® (t) + J [yi (¢— 1) B(t) + y(t — 1) E(t)] dt. 


A material with this symmetry is called hemthedral isotropic. A sugar solution 
containing unequal amounts of right- and left-handed sugar molecules may serve 
as an example of a hemihedral isotropic material. The examples (2.6) and (2.7) 
suffice to indicate how one may proceed to obtain explicit representations for 
the linear functionals (2.3) for the case of an arbitrary symmetry group {S}. 


3. The propagation of an infinite plane electromagnetic wave 
in a hemihedral isotropic material 

An infinite plane progressive vector wave A (a, f) is a time-dependent vector 
field having the form A(x, t) =Bacitm-x—or), (3.1) 
where a=a*+ia~ ({=|/—1, a* and a@~ real vectors) is a complex vector called 
the amplitude of A, k=k*+1k~, k*>0, is the complex wave number, n is a real 
unit vector called the direction of propagation, and w>0 is the angular frequency 
measured in radians per unit time. The symbol # denotes the real part of the 
complex quantity placed after it. 

For each fixed value of #, A is a periodic vector function of time and the 
locus of A is an ellipse or one of its degenerate forms, a circle, or a straight line. 
The major axes of the ellipse, the radius of the circle, or the length and direction 
of the straight line are determined by the complex vector amplitude a. If 
a-a—a’.a*—a.a-+27a*-a-=0, the real and imaginary parts of @ are 
perpendicular to each other and of equal length, and the locus of A is a circle. 
If a*=a*—va@ denotes the complex conjugate of a, and a* Xa=2iat Xa-=0, 
the real and imaginary parts of @ are parallel and the locus of A is a straight 
line. In all other cases, the locus of A is an ellipse. 

In electromagnetic theory (cf. ConpDoN & OpIsHAW 1958, Chap. VI), if the 
electric field E(a, t) is a plane wave, then its amplitude e determines its polari- 
zation. An electromagnetic plane wave is said to be circularly, linearly, or 
elliptically polarized according to whether the locus of E is a circle, a straight 
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line, or an ellipse*. Circularly and elliptically polarized waves are further 
classified as left-handed or right-handed according to the rule 


>0 right-handed, 


at. (a xn) =[a*, a7, n] <0 -— left-handed. 


If a. n=0, the wave is called transverse, and if axn=0, the wave is called 
longitudinal. All other waves will be called skew. 


The field equations which follow from MAXWELL’s equations (1.1) to (1.4) are 


o% + curlE=0, divB=0, 


(3.2) 
curlH—“O —J,, divD=Qp. 


In this section, we shall consider plane electromagnetic waves in a non-conducting 
hemihedral isotropic medium for which J; =Q,;=0, and for which the linear 
functional constitutive relations (2.7) are valid. Thus we seek plane wave | 
solutions, 


E(x, t) = Bee en-e—ot) 
Bt) =F be '" 772, 
D(x,t)=Ad en x—ot) (3.3) 
H(x,t) =Bhe treed), 


of the system of integral and differential equations (2.7) and (3.2). 
Introducing vector fields of the form (3.3) into (2.7) and (3.2), we find that 

they will constitute a solution if and only if the amplitudes e, b, d, and h, the 

wave number k, and the frequency w satisfy jointly the system of equations 


tbe Bini solid cohsista (3.4) 
knxe—wb=0O, (3.5) 
knxh+od=0, (3.6) 
d=a(w) e+ y(o) b, (3.7) 
h=B(o) b+ d(a)e, (3.8) 
where i) Ly 

a(w) = Dinas Oe fF @, (u) eb?" du 

ylo) = Do(— ia)" + foalu) fer du 
pe 2 (3.9) 

B(o) = a b,(— 10)" + J v(m) e°"dy 

d(o) = Ya, iq)? +f yal) Ct ay 


x Linearly polarized electromagnetic waves are also called plane polarized waves, 
but this seems to overwork the word ‘‘plane”’. 
Arch. Rational Mech. Anal., Vol. 6 14 
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From (3.4), (3.7), and (3.8) it follows that E, B, D, and H are transverse. Since 
we assume k*>0, w>0, so that the wave travels in the positive direction of n, 
equations (3.4) will be satisfied as a consequence of (3.5) and (3.6). Thus it is 
sufficient that we consider the four homogeneous vector equations (3.5) to (3.8) 
in the four unknowns e, b,h, and d. In order for there to exist non-trivial 
solutions, the determinant A(k,@) of the 12x12 matrix of coefficients must 
vanish. Thus the complex wave number # is some function of the angular 
frequency w, and this function, whatever it may be, is called the dispersion 


Right circularly polarized B Left-arcularly polarized B 
Fig. 1 


relation. To determine this relation and other properties of the solution we 
proceed as follows. First eliminate d and h from (3.6) using (3.7) and (3.8) to get 


knx(Bb-+ de) +a(ae+yb) =0. (3.10) 
Now eliminate e from (3.10) using (3.5) and m-e=0 to get 
(6 —« $,)"xb r(Z) Ly) b=0. (3.11) 


The scalar product of (3.11) with 6 then yields the condition 


(6+y)b-b=0. (3.12) 


Thus eather the coefficient d0+y=0 or the magnetic flux B must be circularly 
polarized. But if B is circularly polarized so also are E, D, and H, as follows 
from (3.5) to (3.8). 

Let us assume that (y+ 6) =0 so that B is circularly polarized. Since it is 
also transverse, bt, b~, and n are mutually orthogonal, so that one of the two 
geometrical relations in Fig. 1 holds. 


From Figure 1 it is easy to see that 


Pe reel ea b for right-circular polarization, 


+76 for left-circular polarization. 


Using this result and (3.11) we then see that the dispersion relation is 


ww? = w B i B ==; (3.13) 


where the upper sign in the linear term holds for right-circular polarization and 
the lower sign for left-circular polarization. 
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Thus far we have imposed no restrictions on the real constants a,, b,, c,, d 
or on the real memory functions ¢,(u), @2(u), y,(w), and y,(u), other than the 
order condition (2.4). We shall now assume that k*>0 and k=O. It has 
already been pointed out that the first assumption implies that the wave travels 
in the positive direction of m. Subject to k*>0, the second assumption implies 
that a wave of frequency w and small amplitude does not build up, resulting 
in an unstable situation. We may then express the root of (3.13) which satisfies 
k*>0, k =O in the ee 


Vv? 


es 1 a dy BB +6 
= 3B (3.14) 
where the upper sign in the bracket holds for right-circularly polarized B and 
the lower sign for left-circularly polarized B, and the square root taken in each 
case denotes the determination for which k*>0, k=O. 
The wave length is then given by A= 22/k*, the speed of the wave by v=a/h*, 
and the absorption coefficient by k~. From (3.14) we see that, unless the quantity 


9(c) =s (5) (3.15) 
B 

vanishes, the speeds of the right and left-circularly polarized waves will be 

different so that the medium is oftically active (see DRUDE [1902]). Unless the 

quantity 


x (w) =a?) (3.16) 


vanishes, the right and left-circularly polarized waves will have different ab- 
sorption coefficients, as has been observed in some real materials. 


4. Comparison with the Born-Huang theory 
Beginning with a model of an ionic solid consisting of a set of charged mass 
points, the dynamics of which may be treated either classically or quantum 
mechanically, Born & Huane [1954, 1, Chap. VII] derive a relation between 
d and e having the form 


dz-e(k,w)-e, k=hn, (4.1) 
where the dielectric tensor € is Hermitian, 
€:;(k,@) = ef (k,@), (4.2) 
and satisfies also the relations 
€;;(k, @) = e7;(— ko). (4.3) 


Thus if we put 
e=—e +e, 


the conditions (4.2) and (4.3) are satisfied if and only if 


ej; (Ke, @) = gj; (Fe, @), 

€j;(k, @) = — gi (Kk, @), (4.4) 
ej; (k,@) = e;(— ko), 

6; (k, w) = — &;(— Kk, @) 


14* 
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Now from our equations (3.5) and (3.7), we get a relation like (4.1) with 
6; (Kk, 0) = 00 ;; + 2" einjme, (4.5) 


where e; ;, is the completely antisymmetric axial tensor with e¢;23=1. Thus we get 


Ej = a" O55 + eh Cin 7 Me» 
(39) 
ne (4.6) 
Ej; — a 0; a —. = Chg nN; - 
The four conditions (4.4) are satisfied if and only if 
(aA fe) oath Neg eee Oe (4.7) 


- 


It is obvious that these conditions are not met by an arbitrary assignment of 
the constants a,,b,,c,,4,, and the memory functions 91, Y2,%, and pa. In 
the special case of our relations for which «(w), B(w) are real and positive and 
the functions y() and 6(m) are pure imaginary, & is real, there is no absorption, 
and the conditions (4.7) are both satisfied. 


Further comparison of the two theories is not warranted since even the 
partial overlapping of results outlined above is hardly more than accidental. 
In the Born-Huang treatment, absorption is represented by a transfer of electro- 
magnetic energy to mechanical energy of lattice vibrations, while in the simple 
phenomenological theory presented here, any motion or deformation of the 
medium has been ignored; moreover, we have assumed that the medium is 
homogeneous, while the treatment of BorN & HUANG is based on a highly 
inhomogeneous material medium, namely, a system of mass points. 
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1. Introduction 


Under certain operating conditions the stator and rotor blades of axial-flow 
compressors may show the phenomenon of flutter. These oscillations may cause 
the breaking of the blades. One important task in the study of flutter is calcu- 
lation of the unsteady velocity field and of the forces acting on the oscillating 
blades. For the single airfoil this problem has already been solved satisfactorily 
by various scientists, e.g. [2, 2, 3] and others in numerous investigations. Owing 
to the complexity of the problem in the case of cascades of airfoils, various authors 
on this subject have obtained solutions only under strongly simplifying assump- 
tions. 

The annular channel where the stator or rotor blades operate is assumed 
to be of height small in comparison with the radius, so that radial dependence 
is negligible. The plane cascade resulting from development of a blade row 
is assumed to consist of thin, slightly cambered profiles which in general 
are only slightly loaded if there is only a steady flow. It is supposed that the 
fluid is incompressible and nonviscous. Some papers reporting on the investi- 
gations of these problems are the following. 

BILLINGTON [4], SOHNGEN [5] and CHANG & Cuu [6] evaluated the air forces 
acting on a vertical cascade of arbitrary spacing when the profiles make 
small harmonic oscillations in phase. They used different methods. LiLiry [7] 
and LEGENDRE [8] treated the case where adjacent airfoils oscillate harmonically 
180° out of phase. The problem of a single airfoil executing small harmonic 
oscillations and being centered between two parallel walls is equivalent to the 
latter and was solved by Timman [9]. MENDELSON & CARROLL [10] succeeded 
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in obtaining some formulae for lift and moment, if the line profiles of a vertical 
cascade with arbitrary spacing oscillate with a phase shift of 0°, 90° or 180°. 
Woops [11] by using conformal transformation and NickEL [12] by solving the 
appropriate integral equations generalized the investigations, extending them to 
arbitrary, unsteady flows through a cascade of line profiles, the pressure distri- 
butions being equal on adjacent blades. SOHNGEN [13] evaluated the forces on 
a loaded, plane cascade of staggered plates with small spacing, if independently 
of each other the line profiles execute small harmonic oscillations. He succeeded 
in verifying that a loaded blade row — free to move only in bending — can 
flutter. In his thesis Sisto [/4| produced approximate calculations for cascades 
the profiles of which oscillate with an arbitrary phase shift. 

It is of great importance that SOHNGEN & Quick [15, 16] pointed out that 
only in the case of small frequencies and small Mach numbers of the basic flow 
the fluid may be regarded as incompressible. This fact is in full contrast with 
the single airfoil in unbounded space, since there are generalized characteristic 
modes of the gas in the annular channel. Therefore the potential equation 
dominating the incompressible flow problems must be replaced by the wave 
equation; of course, that makes all the calculations more difficult. The forces 
have been evaluated only for a cascade with very small spacing [16]. In the 
future it will be necessary to concentrate all the investigations on unsteady 
subsonic (or supersonic) flows through cascades of profiles. 

Nevertheless, the author believes that the following problem is of some 
interest, because it generalizes those which were cited in [1] to [14]. The problem 
is the following: 

Find the complex perturbation velocity field induced by the harmonic os- 
cillations of the line profiles in a staggered cascade, and develop formulae for 
the unsteady forces acting on the separate blades. The basic flow will be assumed 
parallel to the line profiles in their mean positions. The N blades of the original 
blade row may move independently of each other with small amplitudes of 
deflection. 

We shall solve this problem in the following way: In Chapter 2 we shall 
state the formula for the required field of velocities in the complex z-plane. 
This field will be induced by the distribution of vortices on the line profiles and 
in their wakes. The boundary conditions will result in a system of singular 
integral equations for the vortex densities. According to KELVIN’s theorem and 
the assumption, quite usual in a linear flow theory, that the free vortices remain 
at the place where they arise, we shall be able to evaluate the densities of the 
vortices in the wakes, to within one constant factor each, ¢,, in the case of 
harmonic time-dependence. In the course of our later calculations the e, will 
be taken as known. Then we can evaluate the downwash on the profiles induced 
by the trailing vortices. The system of integral equations for the determination 
of the vortex densities on the profiles can easily be transformed into a new 
‘ system which corresponds to a steady complex velocity field that is irrotational 
in the wakes. 

Unfortunately, direct solution of the boundary-value problem by means of 
the integral equations cannot be obtained in the general case. For this reason, 
in Chapter 3 the z-plane will be mapped conformally onto the ¢-plane in such 
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a way that the staggered cascade turns into a horizontal one. We shall observe 
that the edges of the profiles in the z-plane will not be transformed into the 
edges of the corresponding profiles in the ¢-plane. 

In Chapter 3 we split the transferred boundary values in the ¢-plane into 
an even and an odd part, regarding the two sides of the profiles. The boundary- 
value problems arising may be formulated under the Hilbert type. They will 
be solved by methods described in the book of MUSKHELISHVILI [17]. We shall 
obtain a unique solution if we take into consideration the fact that the per- 
turbation velocities must vanish far upstream and remain bounded at the trail- 
ing edges of the line profiles, due to Kurta’s condition in the original z-plane. 

In order to obtain the unsteady perturbation flow required originally, the 
N as yet unknown amplitudes ¢, of the trailing vortices must be determined 
by the Kelvin circulation theorem. This condition results in a system of linear 
equations for these N unknowns, which will be separated by means of a suitable 
Fourier transform. 

If the perturbation velocity field is known, the lift and moment acting on 
each blade can be stated. In Chapter 5 we shall derive formulae for the Fourier 
transforms of these quantities. We wish to emphasize that these formulae are 
only of theoretical interest; for practical calculations they are too complicated. 
At this point approximate methods must be applied, and a beginning should 
be furnished by the evaluation of the integrals stated in the Appendix. We shall 
not carry the work further, since it would require some lengthy investigations. 


2. Reduction of the Unsteady Flow Problem to a Steady One 


Let us agree upon the notation for the parameters which specify the geometry 
of the cascade. 


a) 2c=chord length of the profiles, 

b) a= distance of the leading edges of the profiles, 
c) I = a/c = spacing of the cascade, 

d) A= angle of stagger, —a/2<AS+a/2. 


In the plane cf the cascade we introduce the dimensional Cartesian co-ordi- 
nates %, Yo and the complex co-ordinate z)=%)+7-¥». and at the same time 
a non-dimensional Cartesian system Oxy. The connexion between the co- 
ordinate systems is given by zy=c-z. The individual line profiles of the cascade 
are given analytically by %=C+ x’ +t-am e~*4 or zo=c-2,,[%'], m an integer, 
and 2,,[%"J=x«'+1-Tme"*, —1<x'<+1. Parallel to its positive x,-axis this 
Staggered cascade of line profiles is approached by an incompressible and 
inviscid fluid at the rate Uj. The individual points 2m[x'] of the profiles 
are assumed to move from their mean positions in harmonic deflections of small 
amplitude. Let the harmonic time-dependence be given by the frequency ». 
In accordance with the linear perturbation theory we assume the profiles to 
remain in their mean positions given above and only the prescribed velocities 
normal to the blades to vary. The deflection velocity of the point z,,[%'] 
on the m'" blade in the direction of the positive y-axis at the time ¢ is denoted 
by Uy-8 q(x’) exp (2707 t) =U,-¢(z,, [x’]) -exp(2avjt), —1<x'<+41. Since 
the cascade is thought of as resulting from the development of a blade row 
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with N blades, 
(4) n(X) =8,(%') for m=k+n-N 
with k=O, 1,..., N—1 and an arbitrary integer. 

Here we have the following problem: Find a complex perturbation velocity 
field with the same harmonic time-dependence, 

U,)- W(z) -exp(2a¥jt) =U)-[U(x, y) —7- V(x, y)]-exp(2a07 0), 
such that the total velocity field t Yor 
Zo =Iytl-Yo=CZ if 


ww (z, 2) = Uy, + Ux Za=rr+i-y : 
; se SSS ~ 
xX W(z)-exp(2z 17 2) ; 


(2 
| / 
Zo = 0°24 [2 vg Uy-G, (x')-e ert 
satisfies the boundary conditions ~142'S47 oe BY A earn 


G) Im, {1 (2, Laas t)} vA 
/ 
= — Up + &m(%') - exp (2207 #) SSS os ae oF 
Yb / he =tae 
on the blades z=z,,[%’], -—1< ms 
x’<+1, m an integer, and the = 
ten GS} x'=tC 1,22 


Kutta condition ae f 
Yo Ym (L')-e pid? Up Em (x')-0 2 7TFE Zo=C-Zyp, L'] 

tw (z, t) =O(1) .or ED Mh Soa 
(4) 7 


as. 22 [1]. 


/ 
In addition the Kelvin theorem, TPA Ld ee -- 
which we shall soon formulate i 
precisely, must be taken into u 
consideration. Fig. 1. Geometry of the Cascade 
(The two imaginary units 7 and 7 must be distinguished strictly. They must 
not be confused, which means in particular that 77=—1 is not valid. In all 


formulae containing the factor exp(2avj7?t) we mean that the real part with 
respect to 7, Re; must be taken.) We shall write the time-dependent factor in 
the form exp(2717?t) =exp(jws) with s=U)-t/c andw=2arc/U). 

Since the prescribed normal velocities are equal on either side of an indivi- 
dual blade, we try to solve our problem by distributing vortices on the blades. 
For the densities of these “profile vortices” we put U)-y,,(x’)-exp(j@s) for 
—14<%'<+1, m any integer, and demand y,,(x’)=y,(x’) for m=k+n-N, 
k=0,1,..., N—1, ” an arbitrary integer, due to equation (1). The total cir- 
culation around one profile is time dependent and will be written in the form 

41 
(5) Uy: ¢- Dy (8) = Uy 6- J yi(2’) dx" exp (ns). 


Any change in the total circulation around the m' profile causes free vortices 
to be shed from the trailing edge z’=c - z,,[1] according to KELVIN’s circulation 
theorem. We assume that the free vortices are carried away with the basic flow 
at the rate U, into the wake, given by 2’=c-z,,[x’] with +13%’<oo, and 
that their strength remains constant. If the harmonic motion of the profiles 
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has lasted for a very long time, the free vortices will have formed a plane vortex 
sheet in the wake of the m‘* blade. Since the terms for perturbation velocities, 
vortex densities, etc. always contain the factor exp(j@ s), we may cancel it in 
all those equations and restrict our investigations to the amplitudes. In order 
to build up our complex velocity field we start with vortices of strength c - Ug x 
6,(x’) dx’, which are placed at the points z’=c-%[x’]+c-17- N-e-*4, » an 
arbitrary integer. The velocity field induced by such a series of vortices is known 


to be [13] 


oi id ° 
(6) dW, (2, x!) = 6,(x') dx’ - 2 ctnh {ay (z — 2 [x i}. 
However, far upstream this field does not vanish, for 
- , Sf, , } vA 
(7) lim dW) (2, «') = — 6, (#') dx Say 


Our cascade being thought of as resulting from a development, the field far 
upstream must vanish. Thus we add the homogenous field 6, (x’) dx’-7 e'7/2T N 
to the field of equation (6). If y,(x’) is the strength of the “‘profile vortices’’, 
and if ¢,(x’) denotes the strength of the free vortices in the wake of the k™ blade, 
we integrate over the chord length and the wake, respectively, sum up the N 
vortex systems, and obtain the total complex velocity field as follows: 


Dania, ty 
W(z) = Sa V(x") = letnh ( = (2 — 2, [x'])) _ 1) aa 
k=0 ° 
(5) fe 
jie 4 ea 
aoe ex (2’) [et h (Fy (@— ale’) +1] ax 
h=0 t 


This field vanishes as x=Re,;z—+— co and remains bounded as x—>+00, if 
N-1 © 


» f e(x’) dx’ is bounded. 


k=0 1 
The Kelvin circulation theorem requires that the relation 


41 
d pie ' 7 
(9) ge Uae: f yan!) ax -exp(j@s)+ Uf + (1) -exp(j@s) =0 
=f 


shall hold for k=0,1,..., M—1, and, because the free vortices flow off with 
the basic stream, we have 


(10) Uo: &(x") exp (j.@ s) = Ug & (1) -exp[ja(s — x’ + 1)] 


for x’21. In the case of harmonic time dependence equations (9) and (10) 
imply that 


. . . au 
(11) €.(2’) = —jo-exp(j@ —jw x’) - fy, (&) dé 
ei 
for 424" and k=O0,4) ANS Bor conciseness we put 


(12) Se TOURS A eure e 
=i 
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Now the densities of the free vortices are known except for the N constants ¢,. 
We insert the terms for ¢,(x’) into (8), transform ctnh--. +1, and obtain 


PN ate ex gees -dx’ 
ee ee ee 
ITN ee | Sees bas 2m et zy[ x") 
(13) ae TN TN 
} Ba 2merrz : 
ne Sadi ex TN id 
aa tx Dadi = EE? i 2 Ag Cy | | 
Tie a TN 


The unknown vortex densities y,(x’) must now be evaluated from the conditions 
Im,{Up > W(%m[%'])} = — Up+ 8m (x) for —1<%/< +14, 


m any integer. To the latter conditions we add the Kutta condition W(z) =O (1) 
as z->z,,[1]. Let us write the kernel in the formula for W(z) in the following way: 


ep wade 
1 eta it TN 1 1 
SiN = ary ae ee Sek: Pe ene Tx! 
(14) IP IN es 2nerg pe Zee 4. EALAl Qni mlx’) —2 i: (z m2); 
P TN P TN 


R (z—2,[x']) being holomorphic in a full neighborhood of z=z,[x’]. Then by 
[17, p. 84 seg.] we see that W(z)=O(1) since z—>z,[1] is possible only if the 
distribution densities of the “‘generalized integrals of the Cauchy type” in (13) 
are continuous in x’—1; 7.e., if 


(15) Hig ( A) seve, Ter? fOr Ah IN Sa. 


The boundary conditions for the normal velocities on the line profiles result 
in a system of singular integral equations for the y,(x’) with CaucHy’s principal 
values in —1<x%’<-+1. We have 


were 7: 
&:(*) =— sae: > fm ) K(x, x/;1—k, a) dx’ — 
é N-1 ee 
(16) Sere ey f oF. K(x, x’; b—k, A) ax’ 
a 1 
fOr WP OT 4 
with 
KG, a Ry) 
é cos A exp AX (x= 4') cos J cos {24 mar le v)sin2+T(I—#)}} 
= ee ft 4 aa , a ee te es y 
sinh en 4% «’) cos A +sine{ Beet Die Ld 


In order to abbreviate let us write 


co 


(17) ly, ,(%; a) = — — perf oie Kl, xsl —k, A) dx 
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for k,/=0,1,...,N—41 and —1S*=+1. Then we see that My, ,(%; A) is 
dependent only on the class 1— k =7(N). Thus we may write fy ,(%; A) =h,(%; A). 
Considering the ¢, as known, we can transform the above system of integral 


equations into the following: 
N-1 


ge) = a ER hy, (%:5 Ay =H (%*) = q (2 [x]) 
k=0 ' 


(18) led 
4 \ , ze Be =, / 
Sabra ie f vl=’) Ble x ;l—k, A)dx 
k=0 4 
for 1=0,1,..., N—1. This system of integral equations is also obtained from 


the problem of a steady flow passing through a staggered cascade of line profiles 
which slightly deflects the basic flow parallel to the x-axis, if the normal velocities 
q(x’) be prescribed on the profiles z9=c - z,,[«’] and without any trailing vortices 
in the wakes z5=C-Z,,[%'], «’21. Thus we have reduced the original unsteady 
problem to a steady one. After evaluating y,(x’), which must remain bounded 
according to equation (15) as x’—>+1, we still have to evaluate e, from the 
system of linear equations that results from equation (12). 

From the kernel K(x, x’; 1—k, A) we split off the singular part for /=hk and 
x=4x', and according to [17, p. 74] we then deduce that as x >+-1, 


O (| log |x —4 for R= 
(19) hia es a) =] 0 Mleg le — 1) 

O (1) for cess J. 
From it results that as x >+1, 
(20) g(x) =O(\log la —1)|) for LOA, = oh 


if we assume g,(x) to be Hoelder-continuous for —1<*x<-+1. The reduced 
and steady velocity field U,- W*(z) we require is then set up as follows: 


Neat oe exp nets 7a, 
a(>) 220 ; aay 
(24) CN ie = UENO er ET LENEANS 
=0 4 eXP aay exp EN 


By equation (15) we require the boundedness of y, (1). After the splitting of 
the kernel according to equation (14) and according to [J7], W*(z) as z>2,[1] 
must behave like O(|log|z—2[1]||) for J=0,1,...,N—1. Once more our 
reduced problem may be precisely formulated as a function-theoretic one: 

. We require a complex velocity field W*(z) that is holomorphic in the z-plane 
with the exception of the slits z=z,,[x%']=x'+iT-m- e* for —1S2x/S+1 and 
m any integer. This field shall satisfy the following conditions: 

Ae WG BAN ee Vee yea) for all z 


2. W*(2) =0/(1) aS Reg = % = =o 
=0(1) as Re 2 == % =o 


3. Im; W*(z) = — q,(x) for z=2,[x] +7-0 

in — 1 t— e Lise Opdh eehad ee 
4. W*(z) = 0 (|log |z — z,[1]]l) as z—>z,[4] 
be W* (2) =O (|z—4,[ Vy nas ee 1] 


with OS Cy 4b. 
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3. Conformal Mapping 


It is very unlikely that the problem defined in the last section can be solved 
directly in the z-plane because the system of integral equations (18) for the 
vortex densities y,(x’) is very complicated. We shall try to approach it in 
another way. The z-plane will be conformally mapped onto a ¢-plane in such 
a way that the staggered cascade of line profiles is transformed into a horizontal 
one. For the velocity field of a flat plate cascade in a fluid steadily passing 
through it, DURAND’s book gives some formulae [18, p. 91] which help to obtain 
the required mapping function z=/(f). In the case of the vertical cascade we 
were also able to state explicitly the inverse function €=g(z) [19], but this 
will not be possible here. 

From v. KARMAN & BurGERs in [18] the formula for the conformal mapping 
of the exterior of the unit circle in the «=a«,+7-«,-plane on to the exterior 
of the cascade of line profiles with the spacing 7, the chord length 2 and the 
angle of stagger A runs as follows: 


iL 
20 


(22) 2(a) = |e-**log eee 


= ehAt log = are 


with 0<x*<1. There is a relation between T, x and A [18]: 


F = cos A-log[(1 — x?) +- {2%-cosA+ (1+ x4 + 2x? cos 2A)4t] + 


D 
) +sinA-tan1[2x-sind-(1+x4+ 2x%%-cos2A)~?]. 


This equation cannot be solved explicitly for x in the case of general A. In all 
practical cases this equation requires numerical calculation which takes into 
eccount the restriction that O0<*<1. A list of some numerical values for x, T, 
and / is found in a paper by N. Scuoxz [20]. 

Now we shall map the exterior of the unit circle in the «-plane onto the 
exterior of a cascade of line profiles in the €=-+7-7-plane with the spacing T, 
the chord length 2, and the angle of stagger 90° by means of the function 


(24) bla) =e | tog EO og ete 


2m A/x— a— x 
The spacing t must be evaluated from a relation similar to equation (23), namely, 


7 
(25) ae eee 
where x is known from (23). 
The functions z=2z(a) in (22) and €=€(«) in (24) effect in parametric form 
the conformal mapping of the exterior of a staggered cascade of line profiles 


onto a horizontal one. After eliminating «, we may write 


#) T- cosa | 
IU 


The square root here is to be determined in such a way that it is positive for 
real ¢=é on the upper side of the slit —1<&<+1, where the principal branch 
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for the logarithm is chosen. For further investigations we need the derivative 


dz Hy phtey LRCOSA AG «BATU veal peieacae EGER 
(27) eee TC) |sin : -(sin cola om +tanA}. 
The image points of the leading and trailing edges of the profiles in the z-plane 
may be obtained from dz/d€=0 as 
C,.=—f,+m-t at the upper bank, 


(28) At Aes 
Cp=lot+m-t at the lower bank 


with ¢,=— - sin (sin A-sin = and m any integer. 


Besides this it can easily be proved that the following points correspond to 
each other: 


eG A ee a -sin A (upper bank) , 
C=—ipwand 4=— — = -sind (lower bank). 
In the special case A=z/2 equation (26) reduces to z=€ since t=T then. If 


Too, x0 by (23), and too by (25). It is obvious that in this limiting 
case the mapping functions (22) and (24) will tend to 


oa sale F(a ha) 
(29b) boo(a) =— F(a — 4), 


so that equation (26) becomes 
(29c) Zo (€) =sinA-6+7-cosd- (2 —4)!, 


where the square root is taken positive for ¢=&>41 and the branch cut is 
pat Cea 1. 


4. Solution of the Transferred Boundary Value Problem 


We now transfer our problem into the ¢-plane in order to determine the 
complex flow field W*(z). If O=F[z(f)] denotes the complex velocity potential 
of the flow field to be ascertained in the ¢-plane, the following terms express 
the well known relations between the velocity fields in the ¢-plane and those 
in the z-plane: 


(30) u(E,m) —4-0(8,m) = 20) = 4 = WH) «f'(). 


Since the field W*(z) is required to have the period i7Ne~'4, the field Q(C) 
must be characterized by the period t- N. From equation (27) we easily infer 
that dz/d€>iTe'* as |Im¢|=|y|—>co. The conditions W*(z 2) =O gas 
Re z=x—>— oo and =O(14) as Rez=x—-+00 then change into Q2(¢€)=o0(1) as 
Im€=—-+oo and =O(1) as Im¢~=yn=—— oo, In transferring the boundary 
conditions Im W*(z)=—gq,(x) for z==z,[x]+7-0 with —1 <*<-+41 and 
1=0,1,..., N—1, we have to take care that these conditions be satisfied on 
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both sides of the line profiles z= z,[x]. On the staggered cascade the prescribed 
boundary values g,(x) are equal on either side of the segments z=z,[x], but 
in general they are no longer so after conformal mapping onto the ¢-plane. The 
boundary values on the upper side of the slits t-/—1<&<1-1+1, 74 =0 must 
be distinguished strictly from those on the lower side of these cuts. On appli- 
cation of (30), we can state the following boundary values: Im {Q(¢) -dl/dz} =: 
—q(x)=—q(%[x]) for 2[xJ=f(é) in —1<x<+4+1 or t-1—-1<&<1-141, 
respectively, with /=0, 1, ..., N—1. From the above boundary values we obtain, 
after some transformations, the following ones for 2(¢)=w(é, 1) —7-v(&, 7%): 


é€ 
sin 2 (=-—1 
v(, +0) =gif(E+4-0)]- S984.) a es Benen 
(sineZ =sint =) 


—sin (2 —i| 
_L cosa ™\y 


: IU 3 IOS 
sin? — —sin2—— 
5 T 


(34) 
v(€, — 0) =4[7( —2-0)] 


j + tan A} =7; (&) 


for t-—1<F <7 + 1-10; 4,.;:; N—1. The square roots:,(...)? must be 
taken as positive. 

We now study the behavior of the newly prescribed boundary values 7; (&) 
on the lines t-1—1<€<t-l+1. If H(z) is any function in the z-plane, which 
need not be analytic nor defined in a region except on arcs, and for which 
the following estimates are valid (m an integer), we have 


O(|z—2,,[—1]]) for z>z,,[— 1] 
O(|z —2,[1]]) for z>2,,[+ 1] 
(32a) H(z) = O(|2—2»[|— sin 4]]} for 2> 2, < sin] +4-0 


O(|z 2, — = sin A]|) for 2 Z,|—— sin dl —4-0; 


then for the function H[z(¢)]| transplanted into the ¢-plane the following estimates 
will hold in the image points: 


O(\C—t-m+C,|*) for €>1t-m—l,+1-0 
O(j¢—t-m—C,|*?) for C>1t-m+l,—71-0 
O(|\C—t-m—1]!) for C>t-m+1 
O(\C—r-m+1|*) for€>t-m—1. 


(32b) H[z(¢)| = 


Since ¢(z[*])=O(1) for x>—1 and =O(|log|*—4]||) as x>+1, /=0,14,..., 
N —1, the new boundary values have the following properties: 


i(E z O(|é—17-2+4)]) for E>1-1t—Zy 
Sore messi for E>1-1#1, 
(33) —_ O(|é —1-1 —6,| - |log | —2- t —Sgl|) fOLae ey eee 

Lee eh eA for F>1-t $1. 
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From the behavior of W*(z) at the points z=z,,[+1], m an integer, and from 
the fact that W*(z)=O(1) as z—>z,,[*] for fixed x when —1<%<+1, and 
from taking into consideration equations (30), (27), and (32), we infer the 
appropriate behavior of the transformed field Q2(¢) as follows: 


O(|€ —t-m—L,|—*e*") for €>t-m—l,+1-0 
00) = O(\C —1-m —Lo| - [log |f — t-m — Coll) for C>1t-m+6,—71-0 
O(|\¢ —t-m +1|-*) for C>t-m+ 1. 
Since OXa,,<1, we write —2a,+1=—M» with OSu,<1, and we know that 


on the slits 7: m—1<é<1-m+1, 7=0, Q(C) has at most integrable singu- 
larities. Since in general the prescribed boundary values 77(€) and 7 (&) on 
both sides of the slits are not equal to each other, we construct the required 
flow field Q() from two others, namely 2, (¢) and 2,(¢), in the following way: 
We define new prescribed boundary values 7,;(&) and 7;2(&) by 


(34) rn(é) =21n(§) +7 (§)] jeri, Wet 


nmalE) = 3 [7 (6) — 7 (4) 


on the upper sides and by 7,,(&) and —7,() on the lower sides of the slits 
t-l—1<é<1-14+1, n=0, respectively. Let the field Q,(¢) induced by vortex 
distributions have the boundary values 7,,(&), and let the field 2,(¢) induced 
by source distributions have the boundary values +7,,(&), respectively. Now 
we specify our problem for the two functions Q, (¢) and Q,(¢): 


We require two complex velocity fields, Q,(€) and Q,(C), that are holomorphic 
and periodic with the period t-N in the C=&+1-n-plane cut along the segments 
Teh Ae 1 Oe cas |7| —>oo these two fields shall remain bounded, 
but 2, (€)+2,(C) shall tend to zero as n->+ 00. When € approaches the real C-axts, 
the following conditions shall hold: 


lim Im 42, (C) les 22, (0) = — 7 (8), 


n> + 
lim ImQ,(¢) =— lim Im,(¢) = — 7,(é) 
n— +0 n—> —0 


for L=0,1,...,N—1. In the limiting points of the segments, C=1t-m—+1, and 
im C=t-m+ly, m any integer, Q,(C) and Q,(f) shall have at most integrable 


singularities. For C->t-m-+y on the lower side it is required further (Kutta 
condition) that 


(2, (0) + £25 (C) | ; : = 0(|log \C.— tm =F Sell) 


It is easily understood that a homogeneous flow field parallel to the real ¢-axis 
can be superposed without any changes in the boundary conditions and the 
Q-conditions on the slits. As |7|->0v, the field Q, (C) +2, (¢) plus this parallel 
field remains bounded. On the other hand, by superposing such a field on the 
field £2, (C)-+,(¢), we can make the total field vanish as n—>-+co, provided 
lim [$2; (¢) +2, (¢)] be real and =-0 as 7—>+00. In the following it is sufficient 
therefore to assume that 2,(¢) and Q,(¢) are bounded as |7|—>co and that 
2; (¢) + 2,(¢) tend to a real limit as y+ 0. 
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We start with the evaluation of the function Q,(¢), since this can be done 
more easily. At the points ¢=o,+”-Nt, n an arbitrary integer and with 
fixed w, in t-k—1<@,<1t-k+1, kR=0,1,...,N—1, we place sources of 
equal strengths Q,(@,) dw,. Then the field of such a series of sources is defined 
by 


dQ, = Q, (a) doy: =, ctn 7 (6 — @). 
For the total source field that results from the integration and summation over 


t-R—1<@,<t-k+1 and over & from 0 to N—4, respectively, we obtain 
the term 


N-1 t-k+1 
(35) D(C) = sta Df Onlos) ct ZC — 0%) deny. 
k=0 7.41 


The still unknown source strengths Q,(@,) must be evaluated from the boundary 
conditions. In a full neighborhood of =a, the kernel can be written in the 
form 


1 cf 1 
Py GA es oe te te oe 


In this neighborhood P(¢—,) is holomorphic. On consideration of this separa- 
tion and the Plemelj formulae [/7, p. 43], we obtain 


(36a) 225(§ + 4-0) — 2,(§ —7-0) =—71- Q(€) 


fort —A< b= 7+l-+41.and /=0,1,.:., N—4. Since 


It 


aN (¢ Wy) , 


ctn a (¢ — w,) =ctn 


it follows that 2,(f¢)=2,(¢). We may then conclude that 


(36b) 22(E +1-0) —2,(§ + 2-0) = —7- Q,(6). 


On the other hand, this is equal to 27 -Im-,(€+7-0)=—272-7.(&) due to 
the boundary conditions. By comparison it follows that Q,(@,)=2-+7,2(@,) in 
t-1—1<a@,<t-l+1. Inserting this result into equation (35), we obtain 


ia tk+1 
(37) Q(Q= Ay df reales) -ctn Sy C — 0) deny. 
k=0 @z=tk—1 


Owing to the behavior of 7 (w,) and consequently to that of 7,(@,) at the limit- 
ing points @,=t-k-+1, the growth properties of 2,(¢) in €=t-m—+1, m any 
integer, in view of [17, p. 74] are as follows: 2, (¢) =O (|¢ —t- m= 1|~#). Since the 
density functions 7; 9 (@;) are otherwise Hoelder-continuous, Q, (¢) remains bounded 


when ¢ approaches the inner points of the slit sides. Since lim ctn a (C—o,)= Ft 
as Imf=7—>-+o9, we have 


(38) lim 0,(0)=F - i PAC Lue 
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Herein we substitute w,=0, (2,[*']) fort -k-—1<@,<t-k+1and —1<4’<+1, 
and after some transformations we may rewrite the last term as 


: N-1 
Jim 20) =F aD fia [9 (2e[x’] + 4-0) —9(e["] —4-0)]dx’=0, 
since the prescribed values Rs rg ted) on both sides of the line profiles 


z=2,[%'] with —1S%’S+1 and k=0,1,..., N—1 are the same. The total 
source strengths in the ¢-plane vanish, resulting from the fact that they vanish 
for the field W*(z) in the z-plane and do not change their values after conformal 
mapping. Using (34) and (31), we insert the terms for 7;,2(@,), substitute 
@,=t-k-+o where te ee and obtain the following terms for 2, (C): 


2, (2) ee fonZe —og—t-k)X 
(39) ‘srt . 106 
sin i 
x | (ge (0) — g (a) «tan A+ (ge (0) + 9% (0) - z|4¢- 


2 a It Oo ' 
sin? — —sin?—— 
E T 


In order to abbreviate we have put gj (o)=q[f(o+t-k-+7-0)],k=0,1,...,.N—1. 
Hence the complex velocity field 2,(¢) is determined. 

Now we shall determine the general solution for the field 2,(¢). As stated 
in [19], we may start with vortex distributions on the slits, and, applying the 
NIcKEL formulae [2/], we may evaluate the vortex densities from a system of 
singular integral equations. In the present case, however, it serves our purpose 
better to reformulate the problem as a boundary-value problem of the Hilbert 
type and apply to it those methods given in [17] that were so useful in reaching 
our goal in [19]. For a field 2,(¢) induced by vortex distributions along the 
slits t-m—1<5§<t-m-+1, n=0, m any integer, we know that 


(40) Q, (6) =—-2,(0). 


Due to the boundary conditions for Q, (¢) we can deduce for t -1—1<&<t-14+1, 
t=O; 1, Sen, NA 


(41 a) 2 (E + 2- = 2, ( —7-0), 
or 2 (& —7-0) + 27-ImQ,(é —7-0), 
or = LOG 4-0) — 24-7, (€); 


between the slits, 7.e., t-m-+1 <&<t(m+1)—1, m any integer, Q,(¢) is holo- 
morphic, so that 


(41 b) Q(E + 4-0) =Q,(— —7-0). 


Equations (41a,b) express an inhomogenous, linear Hilbert boundary-value 
problem for the function Q,(¢). Some additional conditions are 


(41 c) Q(C) =Q2,(€+7-N) forall, 
O C= , ees —Fin ,, 
(44d) ae) =| (| Tm | ") for C>t-m+b1,) OS p,<1, “or 
O(|C —t-m+6,|-*) for €>1t-m—l, 
on the upper side and for +r - m-+-¢, on the lower side of the slit with 02%, ae 
(41 e) $2(0) =O(1) as |n| >. 
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As is known, the general inhomogeneous solution of such a linear Hilbert boundary- 
value problem is obtained by superposition of any arbitrary inhomogeneous 
solution and all the homogeneous solutions [17]. First we determine the totality 
of homogeneous solutions, and for the time being we distinguish between the 
cases of even and odd N. 


a) N 1s an even number. It can easily be verified that the function 


(42) ORC) = Sap aaa [sin C+ 1)-sin ZG —1)] 


satisfies all the conditions made in (41) concerning 7,,(&)=0 and that u,,=4, 
k,,=0. Besides this we obtain 


(43) QC) = O[exp(—F |nl)) as |n|> ov. 
The square root in (42) is to be determined so that 


: : a ee : : mw is 
sin 2 (¢ + 1) ‘sine (fe 1) = 1: isin? Susin2ns. 
1 tT = wv €. 


with the second root being positive for €=&+7-0 and —1<&<+1. The slits 
t-m—15E<1t-m-+1, 47=0, m being any integer, are the branch cuts of rey (¢). 

Assuming X(¢) to be any other solution of the homogeneous Hilbert problem 
with the properties of equations (41), we write the quotient 

_ *<(?) 

(44) HC) = Bey 
Concerning all real points €=é with the possible exception of €=t-m-+1 or 
C=t-m+l,, m any integer, this function H(¢) has the property H(€+7-0)= 
H(é—7-0); 1.e., it must be meromorphic and also periodic with the period 
t-N due to (41c). The conditions (41d) require of H(¢) orders of infinity which 
are less than 1 in the above exceptional points. These, however, can only be 
removable singularities for a meromorphic function. Thus H(¢) is known to 
be an entire function with the period t- N and can be expanded into a Fourier 
series which is convergent in the whole plane: 


Bah exp (22. we) 


v= —OOoO 


Since 2° (@) = —09(0) and X(Q) = — X(0,, it follows that H(¢)=H(¢) 
in a relation among the Fourier coefficients of H(¢): a_,=4d, for esa ft Be Se, 
We required further that 


X(¢) = HC) -Q2(Q)=O(1) as |n| >. 
By means of eee, (43) we obtain 


Say [ \ 
i ye exp (— 2n 2th) <¢ as |n| oo. 
This, however, is possible only if a,=0 for |y| 2N/2+1. 
We see that H(¢) is reduced to a finite Fourier polynomial: 


+N/2 ae i: 
= »' a,-exp (205 =] with a_,=4, 


v= —N/2 
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containing the free complex constants a,. Then the general solution of the 
homogeneous Hilbert boundary-value problem necessarily has the following form: 


X() = RPA [sin 2 +4) -sin 2-1) "x 
(45) 


N/2 


|e0+ D6 cos (202) +d, sin (27 , rel 


with N-+1 arbitrary real constants. It is easily verified that all those properties 
mentioned under equations (41) actually hold for the functions X(¢) with any 
real ¢,, d,, in particular with u,,=%, k,,=0 as here. 
Now we construct an inhomogeneous solution. For this purpose we consider 
Q9(E+i-0) — for T-m—1<&é<t-m+1 
23 (§—7 + 0) +41 for t-m+1<&<t(m+1) — 
This we insert into equations (41a, b), respectively, make a little transformation, 
and then obtain 


(46a) A(E+i-0 Bilt ato [meh 1 ipl pe pee ee, 
1(€+72 : 0) iM 0) (0) in t-2+1<€<t(l+1)— 

Further, we deduce from results on page 210 that 

oe NEE) eater alt ‘~ 


Q(C+ tN) — 28 (C) 
Due to the splitting of the ctn-kernel as described earlier and according to the 
Plemelj formulae (loc. cit.), this inhomogenous Hilbert problem is solved by 
t-k+1 


— 21 + rp, (Wp) Tt 
Q(x +i - 0) ctn Nic (¢ Wz) AW,. 


New 


O,=Tt-k— 
Using (34) and (31), we substitute the terms for 7,,(@,) and consider 

Oo 1, he = (— k+l, ie cae . 
(47) 1(©; + 7-0) = (— 1) Ae 


with the positive square root; next we put w,=o+1-k. Then we are able 
to state the general inhomogenous solution of the linear Hilbert boundary value 
problem with N+-1 free real constants as follows: 


. 7 Ww, Sa 
sine 
t 


. IU 
Sie 
on 


yal a 
xy (—1)* f ctn (Oe ge eRe 
k=0 o=-1 
sine 
(48) x | (94 (o) + qi (0)) - tan A+ (gt (0) — ge (0)) - 
(sin? ae oak sin? ey 
Le HM % 
ah a T-cosh |. -4 
x (sin 7 — sin? 22) do +S sin@ aad +4) - sin ¢—1)| x 


ME 
x \¢ Te, Ch eas (22 He) + d,-sin (2 a2 


a 
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As €—>tm +1 and€->t- m+, this general solution has the required O-behavior 
with the special exponents w,,=%, k,,=0, for the functions g{(o) are Hoelder- 
continuous everywhere and the gj (a) fail to be Hoelder-continuous at just one 
point apiece. These latter have logarithmic singularities at o=¢,, which corre- 
spond to the trailing edges z,[1] of the original cascade. In the formula above 
the g;,(c) are multiplied by the function tan A (sin? = — sin? a2)" sin =? , which 
vanishes at these points in first order. ‘ é ‘ 

Now the remaining conditions for the total field Q(¢) =, (¢)+2,(¢) must 
be satisfied. Let us begin with the requirement Q(¢) 0 as n>+0. Now 


: it's A : 
(49) lim 4 (2) = — 5 (enya + 4 dja) - 


>= CoO 


We substitute the new field Q, (¢) +-Qj (¢) for the old field Q, (2). Since lim Q, (£) =0 
as 7>-+ co and so as to avoid any change of the boundary conditions, we put 


T cosh 
The WN real constants, ¢y,..., Cyjg and d,,...,@yjg_1, are still free; they are 


evaluated from the conditions 
[, (0) + QO) +.:2,(0]-S2 =0(|log |¢ — r-m — ll) 
dz 


as€—>t-m-+,—1-0, m being any integer. By the aid of the Plemelj formulae 
[loc. cit.| and (48), (39), and (27) we have 


oe NAS cee ee) a Gs A ts sin 7 +tan/- (sin? * =—'s1n? =i x 


1 N-1 
Oth 
x (sin? —sin?” aie 
4 if 


ae (sin? = — sin? 2 aa Cn j2— (— 1)! yD (— 1)" fale +r) x 


+1 
SNS fe a(o + TR) ott ago oe o—t-k)do+ 


o=-—l1 


(54) 


x (sin? a sin? ae ctn <a (§—ao—t-k)do+ 


+ (— 1} 60 +e “cos ( ws) + d,-sin (2a + =| + Cyjz* COS al 


= O(|log|& — t- 4 —Zyl)) 


for r-J—1<&é<1-1+1 and €>1-1+€,, J=0,1,..., N—1. %,(o+7- hk) and 
Yp2(6 +t-k) are evaluated by means of (31) and (34). By equations (20) 
and..(32), as €—71-1-+-Co, that, is; as. 2,[¥]—2,[1], the’ first term jis 
O(|log | — + -1—{ |). The factor in brackets has a zero at 6=1-1+, by 
(27) and (28). Therefore the O-condition for the total term can be satisfied 
only if the N real coefficients co, ..., Cyjg and d,,..., dyj2_y satisfy the system of 
linear equations that results from causing the term in the large bracket to 
vanish for =1-/+¢, with /=0,1,...,N—1. This system of equations for 
the c,, d, can be transformed into one with separated unknowns by multiplying 
the J* equation by 1/N-exp[—2ain(t-1+,)/Nt] and then summing over /. 


214 E. MEISTER: 


This is done for »=0, 1, ..., N/2, and 


sey no no) 
Co + >! ¢, + cos (2% Na) + d,,+sin (2x Ne) 

n=1 : 
N-1 at (sine = —sin? ae 
k=0 pe ne sin - (o—C,) -sin Nt (—t:k—o) 

2 [sine (in 24) Cre toe wi gee 

(52) 

sin 1¢ tk) —(N+1)o+N-C,}— 

— cos — (a — fp) sin a ae eT -do + 
N-1 i sin J—>_ (N—1) (¢—1-k—o) 

18. Sy (— Ay f rya(ote-®) -sin ~ cosd tr = do. 
i=0 ae, ss sin (o—£,) «sin wr ((—1t-k—o) 


In the neighborhood of o=¢, all the integrals occurring in the preceeding for- 
mulae are of the Cauchy principal value type. Noting that dyj2=0, we sub- 
stitute these last results into formula (48) for 2,(¢). Finally, the source field 
Q, (¢) of formula (39) and the homogenous parallel field Q}(¢) are added. Putting 
all this together, we obtain the field 


1 oH oe EE =i: 
Oe Vem wie [sin = (¢- 1) -sin = (¢—4)| x 
Nei tt ; 
x“ tilottk (sin? 2 — sin? 22)! x 
2a: i (6+ 7A) (sin? 2 7) 
oe iN 1)C+o+7-k—N -C] cos = 
4 =e — = = =| ——_ ot -. =" 
sin—- (¢—1-k—a) sine eee = 
Ty (6 oO sin = -cos A-+-cos 7 
—7r,.(0-+7-k)-sin™-cosax 
t 
(53) Paes * on 
sin yy [((N—1) (¢—o)+17-h] cos — (o—€,) - cos = 
5 OL et er eee T Tt do 
ees y "SR oeee ate daca 4 
sin— ay (6-e—t h) sin * cos A+-cos ae Ly eee (a—Lp) 
45) 


7; et 
+E" 2. ih tra(o +7 k) X 


=—1 


sin - cos 4+ cos = (o—Cp) 
I : — (o—L, 


we IG, 14 
sin—_* Cos A4+cos —€C, 
T 


(si ot nds ze) 
514 —— —sins ——— 
seit ba Gata ee u ag 


j mC : 4 : 
Sin * co; | y : = 
pose 1-+-cos . sin — (o—{,) 
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We evaluate the transformed prescribed values 7,,(¢+71-) and 7,.(o-+1-k) 
from the g[f(o+7-0+17-k)] by means of the formulas (31), and (34). 
The terms in the big parentheses following the integral sign have a zero term 
of first order in o=€, for all €, so that the zero term of the denominator which 


results from sin = (o —€,) is removable. The above-mentioned formula (53) for 


2(¢) corresponds to formula (86) with 2=0 in [19]. In this paper the total 
circulation is assumed to vanish, so that our 2(¢) must vanish as inl oo, In 


the present paper we see additional terms, vzz those with the factor sin ~ - cos A+ 


COS 22 mo in the denominators. 
1A 


b) N ts an odd number. In the case when N is odd, we do not intend to state 
every single step of the investigation, since our approach would be the same 
as that already shown. There is no change at all in the source solution of the 
field Q,(¢). In order to find the fundamental solution of the homogenous Hilbert 
problem with equal boundary values on both sides of the segments, we again 
choose 2% (¢) by equation (42). Now we must take into account the fact that 
only this is correct: Q2(¢+7-N)=—2}(¢). In order to obtain the special 
solution of the inhomogenous problem, instead of the ctn-kernel one must use 
[sina (€ —t-k—o)/Nt]*+ as the kernel. The general homogenous solution which 
remains bounded as |7|—>oo runs as follows: 


X(6) = Sy fe ged) sin ( — Ae 
(54) (N-+1)/2 foe re 
x De c,- COS 2n(v— 4) | + d)-sin j2z(» — >) al 


with N-++1 arbitrary real constants c,, d,,. These ae can be determined from 
the requirements that 2(¢) +0 as 7 ++ co and Q(¢) - aa” (|log | —1-1—Z,| , 
as €>1-1+f,—7-0, l=0,1,..., N—1. Adding Q, 1) Q,(¢), and Q1(¢), w 


again obtain {2(¢) by (53). Thus we no longer need to distinguish between even 
and odd N. 
The prescribed boundary values g; (a) being known, we have solved our 


problem stated on page 208. Its solution is unique. 


If we let 7 tend to infinity, we may conclude from (23) and (25) that x0 
and t->co in such a way that aH. The mapping function (26) will then 


become that given by equation (29c), which may be solved for ¢ by ¢= sinax 
Zoo tcosdA- (1—22,)* with positive square root for z.,=%o+1-0, —1<%o<+1. 
Multiplying Q(¢) in (53) by d¢/dz and letting 7 tend to infinity, we obtain the 
formula for the complex velocity field with a single line profile oscillating in 
an incompressible fluid flow, also neglecting the field induced by the trailing 
vortices as we did in the evaluation of ay thus 


gwen tine f VS Re 
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5. Determination of the Constants of the Trailing Vortices 
Now we shall solve our original problem in the z-plane. By equation (48) 
the functions q,(x) and consequently gj (¢) still contain the unknown quantities 
&,, k=0,1,..., N—1, which we now evaluate. By use of the Kelvin circulation 
theorem we found the relation 


ot 
E=—fo-d°- ye ax . for l==0,1, Vea, 
ast 


The densities y,(x’) of the profile vortices occurred in the expression for the 
reduced velocity field W*(z) in equation (21). Applying the Plemelj formulas 
(loc. cit.] after separating the kernel in (21) into a singular and a regular term, 
we obtain 


(56) W*(z,[x"] 44-0) — W*(z,[x’] — 7-0) =i (¥) 


for —1<«*/<-+1 and /=0,1,...,N—1. Thus the above formula for e, may 
be rewritten as 


6 oe f [W*(z[x’] +7-0) — W*(z,[x'] —7-0)] dx’. 


The complex velocity field W*(z) is holomorphic outside the line profiles z= z,[x’] 
of the staggered cascade and has only integrable singularities on the profiles. 
By the Cauchy theorem the path of integration may be replaced by a simple 
contour &%, which runs around the / profile only, and no more than once, in 
the clockwise direction. Hence ¢,= —7w-e!®- 6 W*(z) dz, which after conformal 


mapping becomes ¢,=—jw-e!®- $Q(¢) dé. A, the image of Y&, is a contour 
1 


Aj 

which merely runs around the segment t-J—1S&<1-1+1, 7=0, and around 

no other one, once and in the clockwise direction. Relative to the imaginary 

unit 7 the integrals ¢ W*(z) dz are real, since only a field with vortices is in 
L, 


U 
question. In evaluating the integrals ¢ Q(¢) d£, we may confine ourselves to 


Ay 

the field ,(¢), since Q,(¢) is a source field with vanishing total source strengths 
and @1(¢) is a homogenous parallel field. After substituting C=o0+17-l, we 
obtain 


(57) g=—jo-e-¢Q(e+7-l)do for 1=0,1,...,.N—1. 
Ay 


After inserting the term for Q,(0-+-1-/) from (53) into (57), we obtain a system 
of linear equations in the ¢,, which we shall not write down here because it 
is too lengthy. By means of the following transformation we shall succeed in 
separating these equations according to the unknown quantities. This trans- 
formation is based on the following conjugate equations: If 10(%)) 2 dare 
be any N functions of x, these functions can always be written in the form 


N 
fe (%) = > ees) -exp(2zi 77) with 


TA) -exp(— ani He for Fh, w=0,4. 2... Nae 
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Therefore we shall transform the system of eel equations for the e, into 
a new one by multiplying the J equation by + --exp(— 2m ie] and adding 


the resulting equations. The results of some eee occurring herein are 
listed in the Appendix. After application of those formulae, the transformed 
system has the form 


jo N 7+ gO) Nee 
—jo bu 
4 — : = g (sin? 2 — sin? 22)" nef 2 1p1(0 +T-k) xX 
! A, = 
TQ 
vs file sin — (o—¢,) CGS 2% 
(59a) x (sin?— — sin ‘ ! 
5 1 BB 
sin — (9—o) sin — cos A+ cos ——* 
Tt Tt tT 
It 
sin — cos A or (GG) He tate 
— (0 +t-k)-— - COS . 
fn, OE zt T na 9s 
sin - Cos fied sin — (o—Co) 


e179 Nt &™ aby - 9 %0\—% 
i — G (sin — — Sin =e) x< 
Jo 1F i 


dl ed N ; 
‘ ers ea A een Fe Le ; 
x ae exp ani (5 n| No 2% ta | her (G+ 7: R) x 
cop ve 
no\b oe (6—65) OTC dodo 
x (sin? sin? _ —--%,9(6 +7-) sin = -cos A} — ee 
sin — (@—9) sin —— (o—Lp) 
Oe = 4ly Do soon NY == 41, 
Since 
ie ape ’ ba \eaat ee é (—2 -kn 
AE ae (o) exp ( 2764 sl = “2, 8 (o) exp | | = 


aS ; 
N k=0 w=0 
N-1 N-1 
1 1 L n[k—u] 
—>- 1 exp(— 24 Fr): D halo exp (— 2m tk 
b= = 


taking account of the fact that hg, (c) =A (c) for k—w=r+m-N, m any 
integer, we obtain the formula 


(60) Gg (c) = pt) (a) —N.  . p= (c) 


for ~=0,1,...,N—1. After changing the orders of the integrations in equa- 
tions (59b) and substituting g*(c) for 7(c) and r$(c) according to the 
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formulas (31) and (34), which were likewise transformed, we may write 


41 
BAG Cae ty 8 i [g* (oc) + ¢- (a) —N- &. (ht (a) + A- (o))] x 


x tan A+[gt(o)-g- (a) —N- €- (h+ (a) —h (a)) | my x 


x (sin? ——sin® ze i o) +Jo2(9) Zz : E | -do 


é oso pile 14 
sin — (a—£y) (sin —-cos A+ cos—— 
T T v 


= [e+ (c) —g- (6) —N - € (h+ (a) —h- (a) ] - tan A+ 


+ [gt (a) + g-O(6) —N- & (h+ (6) + h-O (a) | : 


sin = -cos 4+ ctn = (o—Co) 
T 
ee 
sin cos A--cos —" 
T 


and 


em) 
in: e-1%-2t : : 4 —1 +(n n 
© en ery aa f [sin 2 —0)| [er lap dtegy © (a) =e 


o=-—1 


== INT 6 Eh (h*(g) aL h-™(o))]- tan) a. [gt ig o)— g Oo \— Nes &™ (h*M(g) — 


WG) ) | ae Z i (sin? Sine 22) Ini (9) -do — 


or PE uno 
(sin? — sin? —— 
i Ta 


~ ft)" fete m9 wom - 


a2 h-™ (0) | »tan A ao [gt™ (o) -L ay, (o) on N ‘ ol” ; (At) (o) —- h-™ (0) | x 
ees! 
Sih 


‘ Tt 
Kee Ef e008 A tha (@).d a tutor ear gee head 


These equations are separated ec cota ne to the unknowns «” and can easily 
be evaluated. The functions /:(c) and f?s(c) occurring in the last formulae 
are given in the Appendix. 
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In the preceeding formulae we determine the functions g*™ (co) by means 
of the prescribed velocities of the blade deflections g,(x)= g (% [x]) as given by 


g(a i Delt o+t-14£4-0)]-exp(—20i +7). 


This is equivalent to the separation of the system modes of the blade row into 
harmonic fundamental modes over the circumference, as may be found in a 
paper by F. Lane [22]. For the functions h=™(o) we obtain the expressions 


==(0) —jwx 
io) —— hy [tare 
(62a) as 
‘ get o+70) ae costa T ety w]sin 2} 
sinh®{-7. (f(o+70)—+*’] cos i} +sine{# —[f—+*’]sin i} 
h-O(G) = — cy . ih G 108. lexp CA [f(o +70) —x’]-sin a} x 
i 
x sinh { AD Pea [f(¢ +70) — x’]-cos i} at 
+exp{— eae ast) [f(o +70) —x’] -sin A} x 
: 2 . , 
(62b) x sinh { ae [/(o +40) — x'] cos Ab] x 


ls 
x ‘sinh? iF [7(o+7-0) — x’] -cos Al Bie 


; ; Oe et saw leds 
! sin? {4 [f(o +7-0) — x']-sin }| -d% 
for'#—=1,2....N =. 


The constants « are evaluated by means of the formulae (61a, b); the values 
of the constants ¢,, k=0,1,..., N—1, of the free vortices can be ascertained 
from equations (59). These values we substitute into formula (53) for the total 
flow field in the ¢-plane. However, it serves our purpose better if we make use 
of the Fourier transforms «”, g+™(o), and h=™(c) and rewrite formula (53) 
in terms of these quantities. Because of the relations 


N-1 By 
= £(”) -ex } = ant 
qi. (0) 2.9 (o) exp(2ai Wr) ems Ont Sao ads, 


the summations over & and » in (53) may be interchanged, and the one over k 
can be carried out by means of the formulae in the Appendix. After multiplying 
Q(6) by d£/dz according to (27), we obtain the required reduced complex velocity 
field W*(z) in parametric form as a sum of WN single fields corresponding to 
the N fundamental system modes of the blade row: 


Z Ne : 
Te = 2 2) 


W*[zQ)]=2 
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with z=2(C)=/(C) given by equation (26) and 


4 ree (eee 2 
CHONG) = (sin - sin = x 


Tt \ 
g AG 
ga sin (o—Lp) cos ——- 
ie (0) 2 2 TO 3 t 4 Tt ae 
x | (9) - [sin sin 2 aE irg ae 
sin — (€—o) sin — - cos A-++ cos —— 
ae iG e T 
cos (o—o) - cos es 
Sa one te do 
: T t 
(63a) — 7 (a) - sin - cos A — = — + 
T 5 5 
sin a oe A-++cos ee - } sin (o—£o) 


+1 sin = -cos A = 
vv , 
Be + i 7) (c) je = (¢ —o)4 =—ctn - (ao —fo)¢ — 


; 4 7 
sin — - cos A+ cos 
T 


CT, et 0 \% 
sin? — —sin FA 
T T 


=, ry) (a) 5 = = do, 
; (sin = . cos A-+cos 0) sin — X yt.) 
\ T © 
ag ee f (E—L,) 
QO) (sin? ~ —sin? = j vr") (a) [sin? * = sin? ) : a 
o=—1 - sin — (€—o) 
Dae IN: 
exp — n ( o) 
(63 b) — rf” (o) sin -cosd | vel : do + 


LG, forty — 12 ele 


In the special case when all blades oscillate in phase, all Fourier transforms of 
order »>0 vanish. 


Thus the reduced problem of perturbation flow has been solved completely, 
at least in parametric form. In order to obtain the field of the unsteady per- 
turbation flow in the z-plane, we must add the velocity field induced by the 
free vortices in the wakes z=z,,[x’] with «’21 and m any integer. Finally, 
we multiply the terms for the amplitudes by exp(jws) and then take the real 
part with respect to 7. We state the field induced by the free vortices, W,, (2) 
as a function of the e omitting the time-dependent factor: 


See Qnetrz 

4 a ca NG 

W,,( ae € fe jax’ i Olen : - , 
A > “ exp ges expe we te 
(64) N-1 se AN 


Wale) => WY" (2) 


n=0 
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with 
SS 2m etrz 
eh) : C2 
6 a WwW.) ae vet aKa T F 
( 5 ) eS (2) T sis a ae Qmerr2z QZ eth x’ op 
Soe Leg i] 
iA 
Ay Ge exp = (2n—N) (z— +’) 
WAN" (2) = Sf enter EN dx! 
DIE i A, 
(65b) #=1 sinh 7. (e— 2] 
5 


for #7 =1,2,...,N—14. 


The first term gives the field in the special case when the blades are in syn- 
chronized oscillations. 


6. Derivation of Formulae for Lift and Moment 
Since the perturbation velocity field U, - W(z) is known, at least in parametric 
form, we are able to evaluate the perturbation pressure field 0, - Ug: f(x, ¥, 2) 
with the density @) of the undisturbed fluid. Due to the fact that the pertur- 
bation velocity field must vanish as *—»— oo, we obtain the pressure field 


(66) p(x, 9) =— U(x, 9) —jo- J Ute, 9) dx! 


by means of the linearized Euler equations. For conciseness, we substitute 
U, (x) for U(x+msin dA, mcosA-+0) and in the same way V,=(x) and p45 (4%). 
The pressure jump 44,,(x)=p;,(x)—~,,(x) on the blade profiles z=z,,[x], 
—1<x*<-+1, m being any integer, is then given by 


(67) A Pm (4) == Ynl2) = 70+ Lyle) 22" 


Ym (x) = U, (x) — U,, (x) being the density of the profile vortices. By integrating 
over the chord length, we obtain the lift (per unit of length in radial direction 
of the blade) %,=40,-U}-2c-A,, with 


apa! x 
(68) An=S [yn(#) +e f ym() dx’ |as. 


The moment acting on the m" blade relative to the leading edge and considered 
positive if tail-heavy, is given by M,,=—40)- Up: 2c?- M,, with 


(60) My =—f (x+1)[yml2) +70: Syme!) x! 


Since the cascade was obtained by the development of a blade row, the following 
relations are the consequences: A,,=A, and M,—M, for m=k(N) and 
k=0,..., N—1. By means of partial integration we may transform the formulae 
for A,, and M,, into the following: 


+1 tel! 
(70) Ay= (1+ jo): fyu(x) dx —jo- fxn (sax, 


and 
+1 +1 


et 
(71) M,=—(1+370) L ya(7) dx —(1—jo) eee) dx+ajo-f x?-y,(x) dx. 
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—jo 


; j 
Due to (12) and the definition of the z,[%], we may substitute — oe -€, for 


Syy(x)dx and z,[*]—iTke-"* for x in (70) and (71). Regarding the 


vortex densities y,(x) as the differences of the perturbation velocities from 
the upper to the lower sides of the profiles, i.e., y,(*)=W*(z[*]+7-0) — 
W*(z,[x%]—17-0), we may transform the integrals in (70) and (71) into contour 
integrals around the k" profile. By use of the Cauchy integral theorem we are also 
able to extend the integrations over any simple contour ¥%, which runs round the 


profile just once in the clockwise direction. Then for k=0, 1, ..., N —1 we deduce 
(2). Ap = (=> + ee -E,—] O° b (2g —iT ke’) -W*(z) dz, 
Ly 
M,= (S22 4 Set). — (1 —jo)- Gi The) WH) dat 
L, 
es 4b jo: Ge iT ey. WH ds 
Lr 


By means of the conformal mapping z=/(¢) of the horizontal cascade of line 
profiles onto the staggered one, we can transform the curves of integration Y, 
into simple contours 4, which encircle only the k segment, 7-k —1 S&S 1-k+1, 
n= 0, once in the clockwise direction. If we then put €=7tk-+ 0 and take account 
of the relation {(tk+0)=iTke—**+ (0), we obtain for R=0,1,..., N—1 


(74) Ag=— (SI Fe-F?) jo Hf) - (ot th) do, 


My= (GO + Se). — (1 jo) Gil) -Q(o+ 7-4) dot 


(75) 
+ Zio P@PQe+t-A)do, 


with 2(¢) from equation (53). 

Immediate calculation of the terms A, and M, is not advisable, but it is 
more suitable to calculate the ‘‘Fourier transforms” A“) and M™ that correspond 
to the N fundamental system modes of the blade row. We then apply the trans- 


formation given by equation (58) to the formulae (74) and (75) and, after a 
series of calculations, obtain 


: : 2 
(Ores, sot ~jo NO | oe - > UO\S 
A) = | ae en 46°) = saan | (~” (o) - (sin? = — sin? =?) Ir (0)) -do— 
a 

o=-—1 

fey 3 eee 4 

sie aye ; sin — - cos A - cos — (o—¢,) 
(76a) 12. | 7 (a) | Jos (2) «sin = (6 —£,) + Jo —* : 
oil sin — - cos A-+_cos LT) 
© Tt. 
— ry (9) - [sin sins =)! te er wu 
Tv 7 ) 


4 Tt : Tt 
sin— - cos A+ cos — ep 
= + cos - iy Jf Sia = (o—C,) 


Flow in Cascades 223 


A) ——(< Oe ates ee (72° (0) (sin? = — sin? acy Sui o) = 


J@ : 
o=-1 
+1 
(76b) — 7% (a) -sin 2 - cosd- Ja(o 0)) a i. [ro m3 (9) do, 
sin— (o—€p) cat 
Tt 
1On Hh 12, ou, IN 4 
M) — i a a Sey eh: 4 1 Lf (9 ) -(sin? = — sin? “1 ) teste 
es 
2 
—(1~jo) -fi(o)}-sin = (6 —f) + Fjo-——_#__| _ 
. sin - cos i+-cos oe 
een (ete sin= - cos A cos = (6—£y) do | 
"2 (0) ple? Joe 6 IU Ut eu 
ae sin * cos A+ COs — Gr sin — (o—€p) 
a 
Sui 
os hee , 
a =f 7) (c) sin= (o b)15 j@-Jo3(a) — (1 —jo) - Jos(o)} a 
o=-—1 
sin - cos A+ cos = (o—Lp) 
St OV ee ~ * 
sin = 60s A-+-cos < Ge 
for r: ey TRON: 
(sine i —sin? 25) re 


+ 1(0) (1 fo) - Be: ee a 


5 5 4 IU 
sin — . cos A-+cos — sin — (o— 
: +cos— fy } sin (oy) 


7 @ 2) 
+i 
ow : (rs () - (sin? = sine ney {4 ie Te (c) —(4—Ja) By (a)\ — 
(7b) —r$(o) -sin 2 - cosa {Pfam fia(o) — (1-7) - Iia(o)}) aie. a 
24 
+4. frQ(o) {Fie Se(e) — (tie) Ita(o)} do 
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In the preceding formulae for A and M the zero term in the denominator 
for o=C,, which is due to sin = (a —€,), 1s removable by the zero term of the 


numerator at this point. Thus these are ordinary integrals in a neighborhood 
OL U ==Gy: 

It stands to reason that the above formulae are not very suitable for practical 
calculations of the forces and moments acting on the line profiles of the staggered 
cascade. For this purpose appropriate approximations should be derived in order 
to reduce the extent of the numerical calculations. In this paper, however, 
we shall not enter into any approximate calculations, since it takes too long 
to derive them. 


7. Appendix 
In this section we shall write down without proof the summation and inte- 
gration formulae needed in the last two sections: 


N-1 P 
4 for ~=0 
78a pls Cy: OX a aks Tbe! | =¢0: 
( ) wD, % 2 Ne | rb) * | 0 “for n>O, 
N—1 N/2—1 


are (E+E) 


t 


Cn {1 for =O 


2 alliage a >O;, 


1 : 
(786) nN 2 2, d,-exp| a0e (x-l- co) | -sin | 52” (r-t+6)| 
75C lO C= 


2 f tor #0, 
Ne 
1 I 2uIn 
Bare et exp| 25 (e-1-+65)| 
zit 10 for Mel 
A for lig eee 
2 


(78e) ‘=0 


x ; Cos < (o —¢,) for »=0 
23 N 
sin— (o—£,) | exp Te {(F n)o Stench for n>0, 


(738i) =A 
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‘ Ss sin WE (N—1) e+o+k:t—l1-1—N- cs Leagopt? 
N 
N gee ahh 
1=0 sin +; -(g—o+l1-t—R- a, 
Tt 
(79a) ee 
; sin — = Cee So) for n=0 
sin — (@—o) sin = (9 Cal exp] — ant (x aig 0) ane for m0, 
IU 
, No sin +75 (N—1) (e—o) +t: k—7-1] Sanit 
a=: sin (o—o+/-t—k-T) 
0) LOGee i — 0 
= Diwey ING Dicitek 
exp|— 43% (2 — mn) (9-0) — oe *| LOT 10 
N-1 nl 
a c 2. (sin 2. COs A COS: me) oa a 
: 1=0 
(79¢) cos ~&. (sin * .cosdA + cos io) for n=O 
_ 15 € T 
0) for +0 
1 N=1 Fe 2 jhe 
a oe Cee k)-e" N 
(794) k=0 ks 
j eos— (¢ — a) for \ fr 0 
sine (Co) exp| ae. (= n| (¢ 0)| for +0. 


In the formulae (76a, b) and (77a, b) the following functions given in integral 
Form occur: (0); fio(o), Jis(o); m=0,17..., N—4, and fi; w=1,2. As our 
path of integration A, we may choose the following one due to aes integral 
theorem and due to the fact that all the integrands of the above mentioned integrals 
have integrable singularities in o=-+1. Starting on the upper side of the slit 
—1<Reo<+1, Img=0, we go from o=—1 to the point e=a—e (e>0), 
then along a semicircle of the radius ¢ around g=o to ga=o~+¢ in the clockwise 
direction, and then along the upper side of the sht as far as a=+1. Then we 
go back on the lower side of the segment as far as g=+1. Now we take our 
path along the semicircle of the radius e around @=o in the clockwise direction 
as far as o=o—e and return to the starting point e=—1 on the lower side. 
Finally we let « tend to zero. If the integrand has a non-integrable singularity 
at 9 =o, the integral must be understood in the sense of a Cauchy principal value. 
After a series of calculations which we omit here, we obtain the formulae (|o|<1) 


ede HE TONGS 
(sine 2 se 
E e 


(80a) 4 (o)= b de=0  (Nicxet [23], p. 58), 


A; sin — — 0) 


0 == j AE ee ae = Ba coe 
(80b) Blo) = (sin? = sin? 22) Cos — a0 =2; 


A, 
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exp|— 4% (3 —”) (@-0)| “sin F (0-50 
() = —— : do 
vb [sine sin? =) sin — (g—o) 
ioe) Dae, [fINE 
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=2 i 0, 
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= 2-exp | 57 o(> —n) 


= -P_yix(cos =] ([24], p. 87). 


Here P(t) is the Legendre function of the first kind. Due to the relation 
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cos 2 do 
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with the real logarithm log, ... . 
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Zur mathematischen Theorie akustischer Wellenfelder 


PETER WERNER 


Vorgelegt von C. MULLER 


Finleitung 

Die vorliegende Untersuchung eréffnet eine Reihe von drei Arbeiten, in denen 
Fragen der mathematischen Theorie stationarer akustischer Wellenfelder in in- 
homogenen Medien behandelt werden. Die akustischen Eigenschaften eines 
fliissigen oder gasf6rmigen Mediums lassen sich durch die Dichteverteilung @, (x) 
im Ruhezustand, die Schallgeschwindigkeitsverteilung c (x) und eine nicht negative 
Dampfungskonstante a, deren Bedeutung in §1 erlautert wird, charakterisieren. 
Die mathematische Beschreibung akustischer Schwingungsvorgange fiihrt auf 

die Gleichung 


(4) eoV({-V |= iN a2? @ a oD 4 OF 


c2 Of cat CROs 


die sich aus den Grundgleichungen der Hydrodynamik durch einen Linearisie- 
rungsprozeB ergibt. Hierbei ist Y(r, t) das Potential einer vorgegebenen kon- 
servativen Kraftdichteverteilung & (x, ¢). Die fiir das Schallfeld charakteristischen 
Gr6Ben sind die Geschwindigkeitsverteilung v (x, #) und die Druckverteilung # (r, ¢). 
Sie berechnen sich aus ® durch 


(2) ee 
und } 
(3) p—pby=— SF a+, 


wobei ~, der Druck im Ruhezustand ist. 


Im folgenden setzen wir voraus, daB & (xz, £) zeitharmonisch mit der Frequenz w 
ist, und beschranken uns auf die Diskussion derjenigen Losungen ® von (1), 
die die Gestalt 


(4) P(r, t) = U, (x) coswt-+ U,(z) sinwt 


besitzen. Sie charakterisieren die eingeschwungenen ,,stationdren’‘ akustischen 
Schwingungsvorgange. Setzen wir U=U,+7U,, so laBt sich (4) auch in der 
Form 


(5) P(r, t) = Re{U(z) e 7°} 


schreiben. Die Lésungen ® von (1) der Gestalt (5) entsprechen eindeutig um- 
kehrbar den Lésungen der zeitunabhangigen Gleichung 


: ar(svojee en 
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mit 

(7) 8 = (wm +74) 

und 

(8) We, t) = — — Refi f(x) e**}. 
Hierbei kann x(r) so gewahlt werden, daB 

(9) OS arex(t) <4 

gilt. 


Im Fall eines homogenen Mediums, auf das keine auBeren Krafte wirken, 
verschwindet f, wahrend o, und c konstant sind, so daB (6) in die Helmholtzsche 
Schwingungsgleichung 
(10) AU +22U=0 


mit konstantem x iibergeht. SOMMERFELD erkannte, daB genau diejenigen 
Lésungen von (10) physikalisch realisierbaren Schwingungsvorgangen ent- 
sprechen, die fiir y—co und |x)|=1 die asymptotischen Relationen 


(11) g 


oY 


U(r a9) — 1% Ur me) = 0(-), 


1 

(12) U(r %») =0(-) 
gleichmaBig fiir alle Richtungen ry erfiillen (Sommerfeldsche Ausstrahlungs- 
bedingung) [15]. Wir setzen voraus, da alle betrachteten Medien im AuBeren 
einer gentigend groBen Kugel |r| =R homogen sind und daf& f fiir |x| >R ver- 
schwindet. Dann kénnen die Bedingungen (41), (12) zur Charakterisierung der 
physikalisch realisierbaren Loésungen von (6) tthernommen werden. 

Die Berechnung des von einer zeitharmonischen Kraftdichteverteilung in 
einem unbegrenzten Medium erzeugten stationaren Wellenfelds fiihrt auf folgende 
Aufgabe: 


(A) Es ist eine Funktion U zu bestimmen, die im gesamten Raum der Glei- 
chung (6) und fiir y->oo den Bedingungen (41), (42) geniigt. 


Im Fall, da8 @) und x konstant sind und / hdélderstetig ist, besitzt Auf- 
gabe (A) genau eine Lésung, die durch 


(13) UG) =a ! i’) eitit-8l gy 


dargestellt wird. 


AuBer durch zeitharmonische auBere Krafte kann ein stationares akustisches 
Wellenfeld in einem fliissigen oder gasférmigen Medium M durch einen in M 
eingebetteten festen Kérper K erzeugt werden, der elastische Schwingungen mit 
der Frequenz w vollfiihrt. Auf der Grenzfliche F zwischen K und M stimmen 
die Normalgeschwindigkeiten der Randpunkte von K und M iiberein. Die 
Berechnung des von K und einer zusatzlich gegebenen zeitharmonischen Kraft- 


dichteverteilung / erzeugten stationdren Wellenfelds verlangt die Lésung folgender 
Aufgabe: 
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(B) Es ist eine Funktion U zu bestimmen, die im Au8eren von F der Glei- 
chung (6) und fiir 7+ oo den Bedingungen (11), (12) geniigt und deren Normalab- 
leitung auf  vorgegebene Werte annimmt (Neumannsches AuBenraumproblem). 


Mathematisch gleichwertig, jedoch physikalisch von geringerer Bedeutung ist 
das Dirichletsche AuBenraumproblem, bei dem auf F die Werte von U vorge- 
geben sind. Ihm entspricht physikalisch die Vorgabe der Druckverteilung auf F. 
Die formulierten Au8enraumprobleme sind fiir die Schwingungsgleichung (10) 
von CL. MULLER [//] und R. Leis [7] mit Hilfe der Kapazitatsmethode geldst 
worden. Die in den zitierten Arbeiten dargestellten Ergebnisse erméglichen es 
im Fall eines homogenen Mediums, die durch pulsierende feste Kérper erzeugten 
stationdren akustischen Schwingungsvorgénge mathematisch zu beherrschen. 

Die vorliegende Arbeit und die beiden geplanten weiteren Arbeiten setzen 
es sich zum Ziel, die Methoden zur Untersuchung stationarer akustischer Wellen- 
felder auf inhomogene Medien auszudehnen. Hierzu sollen zundchst die Auf- 
gaben (A) und (B) fiir hélderstetige Funktionen x und / und hélderstetig diffe- 
renzierbare Funktionen @, unter der Voraussetzung gelést werden, daB x und 0, 
fiir |x| >R konstant sind und f fiir |r| > R verschwindet. 

Daritiber hinaus wollen wir Methoden entwickeln, die es erméglichen, die an 
den Grenzflachen verschiedener Medien auftretenden Reflexions- und Brechungs- 
erscheinungen zu behandeln. Als charakteristisches Beispiel betrachten wir eine 
zweimal stetig differenzierbare, regulare, geschlossene Flache Ff’, die den Raum 
in das AuBengebiet G, und das Innengebiet G; zerlegt, und nehmen an, da8 G, 
durch ein Medium M, mit der Schallgeschwindigkeit c,(r), der Dichte ),(r) und 
der Dampfungskonstanten a, und G; durch ein Medium MW, mit den entsprechenden 
Gr6Ben c;(z), 09;(x) und a; ausgefillt wird. Wir setzen voraus, daB die Funk- 
tionen 09,, Q9; und c,, ¢;in ihren Definitionsbereichen G,+ Ff’ und G;+ fF’ hélder- 
stetig differenzierbar bzw. hdélderstetig sind und daB oy, und c¢, fiir |r| >R 
konstant sind. Die Berechnung des durch die Kraftdichteverteilung / in den 
Medien M, und M, erzeugten stationdren Wellenfelds fiihrt zu folgender Aufgabe: 


(C) Es ist eine Funktion U mit folgenden Eigenschaften zu bestimmen: 
a) U erfiillt in G, und in G, die Gleichung (6) mit 


@(@+tda) ga "EG 


oa (t) fiir LEG, : Cq(t)? 

£) \ eens es wi q se } 
oe te iy a cout wien ta)) terilZeG,. 
c;(x)* ; 


b) U geniigt fiir 7—> co den Bedingungen (11), (12). 
c) U ist in G,+F’ und in G;+F’ stetig differenzierbar, und es gilt fir yc’ 


2 
oF 


: ae 3% ut ‘ 
(@,—17) Ula seal le 1) Crt 


ae! eee ou; 
Coa On Qoi On i 


Die Bedingung c), die in § 2 ausfiihrlich erlautert wird, ergibt sich aus der physi- 
kalischen Forderung, daB die Druckverteilung und die Normalkomponente der 
Geschwindigkeit bei dem Durchgang durch die Grenzflache F” stetig sind. Ent- 
sprechend kann Aufgabe (B) verallgemeinert werden. 
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Die vorliegende Arbeit beginnt mit einer ausfiihrlichen Darstellung der physi- 
kalischen Grundlagen, die neben der Herleitung der Gleichung (6) und der zu- 
satzlichen Bedingungen fiir U an Rand- und Grenzflachen eine eingehende 
Diskussion der Energieverhaltnisse enthalt. In § 3 werden einige Lésungen der 
Schwingungsgleichung (10) diskutiert. Insbesondere werden diejenigen Lésungen 
von (10) charakterisiert, die den Ausstrahlungsbedingungen (11), (12) gentigen. 
In § 4 wird gezeigt, daB die Aufgaben (A), (B) und (C) héchstens eine Lésung 
besitzen. Die Argumentation, die in den Eindeutigkeitsbeweisen verwendet wird, 
steht in engem Zusammenhang mit den in §4 und § 2 durchgefiihrten energeti- 
schen Betrachtungen. Die Arbeit schlieBt mit dem Existenznachweis fiir die 
Lésung der Aufgabe (A). 

In der geplanten zweiten Arbeit werden nach dem Hilbertschen Ansatz der 
Parametrix die Neumannsche und die Dirichletsche AuBenraumaufgabe fiir die 
Gleichung (6) gelést, wahrend im Mittelpunkt der dritten Arbeit der Existenz- 
beweis fiir die Aufgabe (C) steht. Die Existenzuntersuchungen kénnen mit Hilfe 
potentialtheoretischer Methoden auf die Diskussion von Integralgleichungen 
zuriickgefiihrt werden. Die auftretenden Integralgleichungen und Integral- 
gleichungssysteme lassen sich unter dem einheitlichen Gesichtspunkt der Theorie 
der linearen Transformationen in Banach-Raumen behandeln. 

Es wird eine moglichst in sich geschlossene Darstellung angestrebt, die vor- 
nehmlich in der ersten Arbeit die Beziehungen zwischen den physikalischen 
Grundlagen und ihrer mathematischen Formulierung betont. Es erscheint daher 
zweckmaBig, auch einige in der Literatur bereits bekannte Dinge aufzunehmen, 
die jedoch bisher noch keine zusammenhangende und unseren Zwecken vollig 
entsprechende Darstellung gefunden haben. 

Fiir die Anregung zu dieser Arbeit und fiir zahlreiche wertvolle Ratschlage 
mdéchte ich meinem verehrten Lehrer, Herrn Professor Dr. Claus MULLER, 
herzlich danken. 


§ 1. Physikalische Grundlagen 


Wir betrachten zunachst ein unbegrenztes fliissiges oder gasformiges Kon- 
tinuum, in dem eine Strémung mit der Geschwindigkeitsverteilung v (zr, ¢), der 
Druckverteilung #(z,¢) und der Dichteverteilung (r,¢) herrscht; v, » und a 
seien stetig differenzierbar. 

Wir setzen voraus, daB die Strémung im Endlichen quellfrei ist. Fiir jedes 
endliche regulare* Gebiet G mit der Randflache F gilt also 
(1.1) — 4, [ ee0aV = f oe.) v(e,y nQyaF, 

G F 
wobei n(x) der ins AuBere von F weisende Normalenvektor von F im Flachen- 
punkt y ist. Nach dem Satz von Gauss ergibt sich hieraus die Kontinuitats- 
gleichung 


a 
(1.2) 5 +V (ev) =0. 


Um weitere Beziehungen zwischen v, p und o@ Zu gewinnen, betrachten wir den 
Flissigkeitspunkt P, der sich zur Zeit t in x befindet. Die Bahn 4(t) von P 


* , regular“ im Sinn von [6], $.112, 113 
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erhalten wir, indem wir die zu den Anfangsbedingungen 


(1.3) a(t) =k 
gehérende Lésung des Differentialgleichungssystems 
(1.4) 8'(t) = (3(z), 7) 


bestimmen. Es sei A(3, t) eine stetig differenzierbare Funktion. Der zeitliche 
Verlauf von A langs der Bahn 3(t) des Fliissigkeitspunktes P wird durch die 
Funktion 


(1.5) F(t) = A (3(2), 2) 


beschrieben. Wir bezeichnen daher f(t) als die substantielle Ableitung der 
Funktion A im Punkt zy zur Zeit ¢. Nach (1.3) und (4.4) ist 


' ) 
(1.6) F(t) =2 AQ,d) +0, RAC.A). 
Fihren wir fiir die substantielle Ableitung das Symbol D/Dt ein, so gilt nach (1.6) 
D 0 
(1.7) ie Spe) 


Fiir die Beschleunigung des zur Zeit ¢ in y befindlichen Fliissigkeitspunktes 


erhalten wir z.B. 
D 


(1.8) ie) Se Mca 
Zur Aufstellung der Bewegungsgleichungen betrachten wir ein Fliissigkeits- 
teilchen w, das zur Zeit ¢ das regulare Gebiet G mit der Randflache F ausfiillt. 
Die in der Fliissigkeit auftretenden inneren Krafte setzen sich aus Druck- und 
Reibungskraften zusammen. Die Druckkrafte lassen sich mit Hilfe der Druck- 
verteilung f(r, ¢) beschreiben, und zwar wirkt auf mw zur Zeit ¢ die Druckkraft 


(1.9) ae p(t, t) u(t) dF = tf V, p(t, t) aV. 


Bei akustischen Vorgangen treten nur kleine Geschwindigkeiten und Dichte- 
schwankungen auf. Wir konnen daher von der Annahme ausgehen, daB auf 
jedes Fliissigkeitselement eine entgegengesetzt zu seiner Bahn gerichtete Rei- 
bungskraft wirkt, die proportional zu seiner Masse und seiner Geschwindigkeit 
ist. Auf « wirkt demnach die Reibungskraft 


(1.10) 2a = a Jolt t) v(x, t) aV, 


+ (vV)». 


™~ 


wobei a eine nur von dem betreffenden Medium abhangige positive Zahl ist. 


Bei der Beschreibung der 4uBeren Krafte nehmen wir an, daB es eine stetige 
Kraftdichteverteilung ®(z,¢) gibt, mit deren Hilfe sich die zur Zeit? auf 
wirkende auBere Kraft in der Form 


(1.41) J (x,t) dV 
G 
darstellen 14Bt. Aus den Newtonschen Grundgleichungen folgt 
2) fel Bd) vee bekag foLP r,t) —Vp(t,t) —ae(t,) v(e,)] av. 
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Diese Beziehung gilt fiir jedes endliche regulare Gebiet G und jede Zeit ¢. Aus 
der Stetigkeit des Integranden folgt daher 


(1.13) opr tagv+Vp=8 
bzw. nach (1.8) 
(1.14) 0 % 4 (wV)o+ar| +7 p=®8. 


Mit der Kontinuitatsgleichung (1.2) und den Gleichungen (1.14), die im reibungs- 
losen Fall a—0 in die Eulerschen Gleichungen iibergehen, haben wir vier Be- 
ziehungen zwischen v, # und g gewonnen. Um eine weitere Bezichung zu erhalten, 
wollen wir die Energieverhaltnisse diskutieren. 

Die Energiedichte setzt sich aus dem kinetischen Energieanteil $@ »? und 
dem inneren Energieanteil 9 e zusammen. Die zur Zeit tin G enthaltene Gesamt- 
energie betragt also 


(1.45) S(ee+ zen?) dV. 

G 
Mit der Flissigkeitsstromung ist ein Energietransport der Starke 
(1.16) f(oe+4ov®) vndF 

F 


durch die Randflache F von G verbunden. Die Energieproduktion in G betragt 
also zur Zeit ¢ 


(1.17) az f (ee+fou)av+ f (eet Sev*jvnar. 
G F 


Setzen wir voraus, daB der Fliissigkeit keine Warme von auBen zugefiihrt wird, 
und sehen wir von den Einfliissen der Warmeleitung ab, so ist nach dem Energie- 
prinzip die Energieproduktion in G gleich der Arbeit, die zur Zeit ¢ durch die 
auf « wirkende duBere Kraft (1.11) und die Druckkraft (1.9) geleistet wird. Die 
Leistungen der Krafte (1.9) und (4.11) erhalt man, indem man jeweils den 
Integranden skalar mit » multipliziert. Somit gilt 


(1.18) fleet Zorjav + [(ee+5ov)ondr— [ xvav—f pondF. 
G F e F 


Nach dem Satz von Gauss folgt hieraus 


(1.19) v —V (pv) =F \e (e+ 2" "a Vor (e+ 4 v*)]. 
Nach (1.2) und (1.7) gilt 
(1.20) vk —V (pv) St (e+). 


Nach (1.13) ist 


D Ne. Di 
e-n7 (5 v*) ov = =v — vp — aor’. 


Hieraus und aus (1.20) folgt 

De 
Wel = =0— — & 
(1.21) PV OO 
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<< 


Nach (1.2) gilt 


ee 1) Oe, es an 
ie e & ree Vo) 
Damit geht (4.21) in 
; De Dies 
> . park paths Se 2 
(4.22) ee ae 


liber. Es sei s(x, ¢) die auf die Masseneinheit bezogene spezifische Entropie im 
Punkt ry zur Zeit 7. Nach dem zweiten Hauptsatz der Thermodynamik kann 
die spezifische innere Energie e(r,¢) als Funktion e(s, 1/e) von s und 1/0 dar- 
gestellt werden, und es gilt 


de oe 


(1.23) ae ae ‘ =—>p, 
= 
Q 
wobei T die absolute Temperatur ist. Aus (1.23) folgt 
IDO ADs ID 4 
(1-24) Dt ah Dt one 


Hieraus und aus (1.22) ergibt sich die Beziehung 


Ds av 
(1.25) >= ry ce 
Die Kontinuitatsgleichung (1.2), die Eulerschen Gleichungen (1.14) sowie die 
Energiegleichung (1.19) bzw. die zu ihr im wesentlichen aquivalente Entropie- 
gleichung (1.25) ergeben zusammen fiinf partielle Differentialgleichungen fiir die 
Komponenten von »v und zwei der thermodynamischen GréBen 4, @, e, T, s. 
Die weiteren thermodynamischen Gr6éBen lassen sich mit Hilfe der Zustands- 
gleichung des betreffenden Mediums ermitteln. 

Bei der Anwendung auf akustische Fragestellungen sind vor allem diejenigen 
Lésungen des Differentialgleichungssystems (1.2), (1.14), (1.25) von Bedeutung, 
fiir die die Geschwindigkeit » und die Abweichungen £,=/— fy und 0, =o — O 
von der Druck- und der Dichteverteilung im kraftefreien Ruhezustand klein 
sind. Hierdurch wird folgende Linearisierung nahegelegt [3]: 

a) Die Zustandsgleichung =#/(o, s) des betreffenden Mediums kann durch 
die lineare Gleichung 


(1.26) Pre = «(0 100) a thems) 
mit 
a. OP eee 
a re) s=const 4 OS |o=const 
O=0Q0,5=So O=Q0, S=So 


ersetzt werden. Hierbei bedeutet s)(r) die spezifische Entropie im Ruhezustand. 
Da der Druck f, im kraftefreien Gleichgewichtszustand konstant ist, folgt aus 


Vp= Vet Ps fir 60, Und s=Sp 


(1.27) ceVo,+AVs=0. 


* Die zu (1.25) fithrende Argumentation ist fiir a = 0 in der Theorie der kompressib- 
len Stromung gebrauchlich. Vgl. etwa [8], S.188ff. 


238 PETER WERNER: 


b) Die Entropiegleichung 


d av? 
(1.28) = +Vs=F- 
wird durch die lineare Gleichung 
0 
(1.29) 5; $V 5 =0 


ersetzt. Eliminieren wir s mit Hilfe von (1.26), so geht (1.29) in 
op ée@ = 


tiber. Nach (1.27) folgt hieraus 


Op _ ,2/e@ 
(1.34) it o/ E+ vV 00). 
c) Die Linearisierung der Kontinuitatsgleichung (1.2) ergibt 
(1232) 22 + 0M 09+ op V0=0. 
Nach (4.31) ist daher 
(1.33) <P + o9c2V0 =0. 


d) Aus den Eulerschen Gleichungen (1.14) folgt durch Linearisierung 
(1.34) Qo-gy + 400+ Vp =. 


Mit (1.33) und (1.34) haben wir vier lineare Differentialgleichungen zur 
Bestimmung von » und # gewonnen. Zu ihrer weiteren Untersuchung setzen 
wir voraus, daB die Kraftdichteverteilung % ein Potential Y besitzt. Es sei also 


(1.35) DeVere 
Aus (1.34) und (1.35) folgt 
(1:36) 3 PX (e08)] + 47x (e92) = 0. 


Es gibt also ein nicht von der Zeit abhangendes Vektorfeld ¢(x) mit 
(1:37) V(65:0) Cee 7% 


Aus (1.37) ersieht man: Verschwindet Vx (gy) v) zu einer Zeit fj im Punkt x, so 
gilt fiir jede Zeit ¢ im Punkt ry Vx(0)9v)=0. In der Theorie der stationdren 
akustischen Wellenfelder beschrankt man sich auf die Untersuchung derjenigen 


Lésungen b, die in jedem Punkt y zeitharmonisch sind. Wir kénnen daher an- 
nehmen, daB fiir alle xy und ¢ 


(1.38) V X(@90) =0 
gilt. Wegen (1.38) gibt es eine Funktion ®(r, ¢) mit 


(1.39) ood =V, ®. 
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Wir zeigen, daB ® so gewahlt werden kann, daB 


(1.40) = -(b-b) —a 


; 5 
gilt. Hierzu betrachten wir mit festem x) die spezielle Lésung ®)= f 0, wtds 


Lo 
von (1.39). Mit t ist auch ®, stetig nach ¢ differenzierbar. Aus (1.34) und (1.35) 
folgt 


(1.41) VS) —V[P — (6 — b) — 4%). 
Es gibt daher eine nur von ¢ abhangige Funktion /(¢) mit 
(1.42) oe =P — (p — py) —a By +f. 
Mit @®, ist auch 

(1.43) P=, + c(t 


eine Lésung von (1.39), wobei c(t) eine beliebige zeitabhangige Funktion ist. 
Ist c(t) speziell eine Lésung der Differentialgleichung 
(1.44) c(t) =— ac(t) +f (2), 
so folgt aus (1.42), daB (41.40) durch ®)+-c(é) erfiillt wird. @ ist durch (1.39) 
und (1.40) bis auf ein additives Glied der Form ce~*’ eindeutig bestimmt. 

Aus (1.33), (4.39) und (1.40) folgt, daB ® der Gleichung 


3 1 1 826 a a 1 OW 
(1.45) QV (V0) ——. or | BOF = ch we 
geniigt. Ist umgekehrt ® eine Lésung von (1.45), so erfiillen 
(1.46) op 

Qo 
und 
(1.47) p—py=¥—-“F —a@ 


die Gleichungen (1.33) und (1.34) mit R=V YW. Fiir den Fall, daB die Kraft- 
dichteverteilung ® ein Potential Y besitzt, ist damit die Untersuchung akustischer 
Wellenfelder im Rahmen der vorstehend entwickelten linearisierten Theorie auf 
die Diskussion der linearen partiellen Differentialgleichung (1.45) zuriickgefiihrt. 
Die fiir die Schallverteilung charakteristischen Gréfen » und # lassen sich aus 
den Lésungen ® von (1.45) mit Hilfe von (1.46) und (1.47) berechnen. Zwischen 
p, s und @ besteht in der linearisierten Theorie der Zusammenhang (1.26). 
Speziell gilt fiir s=sy 


(1.48) P (0, So) — Po = ©? (0 — @o)- 
Ist @9 vom Ort abhangig, so hangt auch die durch 


Ph 
M49) 3 00 s=const 
O=Q0, S=So 
definierte GréBe c vom Ort ab. Im Fall eines homogenen Mediums kann ¢ als 


Schallgeschwindigkeit interpretiert werden, wie in §3 noch ausgefiihrt wird. 
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Wir lassen jetzt die Voraussetzung fallen, daB das betrachtete Medium unbe- 
grenzt ist, und nehmen an, daB Grenzflachen vorhanden sind, die das gasformige 
oder fliissige Medium M, mit der mittleren Dichte @9;, der Schallgeschwindig- 
keit c, und der Reibungskonstanten a, von einem anderen Medium M, mit den 
entsprechenden GréBen Q92, C2, 4, oder von einem festen Korper K trennen. 
An der Grenzflache F(t) zwischen M, und M, miissen die Drucke und die Normal- 
komponenten der Geschwindigkeit in beiden Medien gleich sein. Im Anwendungs- 
bereich der linearisierten Theorie vollfiihren die einzelnen Fliissigkeitspunkte 
kleine Schwingungen um eine feste Ruhelage. Es ist daher méglich, die momen- 
tane Grenzfliche F(t) durch die Grenzflache F zwischen M, und M, im Ruhe- 
zustand zu ersetzen. Damit erhalten wir nach (1.46) und (1.47) fiir rcF die 


Grenzbedingungen 
ec eS 


tS) on on Qo2 én 
und 
PD 
(1.54) oot +a B,— = St + a, By — YY. 


Hierbei bedeutet, falls wir den in das Medium M, weisenden Normalenvektor 
von F im Punkt y mit n(x) bezeichnen, 


Py (t,t) = lim P(x + An(z),?) 


®,(r,1) = lim P(x —An(z),?) 


(1.52) 3 3 
0 Beet 0 : : 
een bs t) Jim, an(e) D(x An(z), t). 


Entsprechend kénnen Grenzbedingungen an den Grenzflachen F’(t) zwischen M, 
und K aufgestellt werden. Wir betrachten den Fall, daB der Kérper K elastische 
Schwingungen um eine feste Ruhelage vollfiihrt. Dann wird dem Medium M, 
an der jeweiligen Grenzflache F’(t) die Normalgeschwindigkeit v, (r, t) des schwin- 
genden Kérpers K aufgezwungen. Nehmen wir an, daB die maximale Aus- 
lenkung der Randpunkte von K aus ihrer Ruhelage klein ist, so kann F’(¢) durch 
die Randflache F’ von K im ruhenden Zustand ersetzt werden. Damit erhalten 
wir nach (1.46) fiir y¢F’ die Randbedingung 


(1.53) 7 P(t, t) = g(t, 4), 


wobei g(r, t) eine gegebene auf F’ erklarte Funktion ist. 

Die auf die Masseneinheit bezogene Energiedichte eines akustischen Wellen- 
feldes betragt 
(1.54) e+ tb?. 


Setzen wir 1/o@=v, so berechnet sich e nach (1.23) durch 
1 


(r,t) s(t) 


(155) e@N=— f P(r, sole) dot f TTL, s\ds-+e(t), 


ae. So (t) 
Qo (¥) 
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wobei é (x) die Verteilung der inneren Energie im Ruhezustand und p= (u, s) 
die Zustandsgleichung des betreffenden Mediums ist. Die innere Energiedichte e 
kann in die Anteile 


1 
s oe 


o— Pol sa | {7 , s}ds und — fp p)av 


Qo 


aufgespalten werden. In der Akustik interessiert vor allem der zweite Anteil, 


der nach RAYLEIGH [12] als potentielle Energiedichte ¢,., angesehen werden kann. 
Wir setzen also 


Q(t, t) 
(1.56) Epot (t,t) = — £ [P (u, So (t)) — po] dv. 
Qo (X) 
Die potentielle Energie 
i Q Epot dV 


eines Fliissigkeitsteilchen w, das zur Zeit t das Gebiet G ausfiillt, kann als die 
Arbeit interpretiert werden, die der Uberdruck £— 4, leistet, um aus dem 
Ruhezustand in den herrschenden physikalischen Zustand zu wberfiihren. 


Wir wollen ¢,,, im Rahmen der linearisierten Theorie durch ® ausdriicken. 
Hierbei k6nnen wir annehmen, daB zwischen # und @ die lineare Beziehung 


(1.48) besteht. Setzen wir me (01 — 09) =41, So gilt fiir 6,0 
0 


—fo—tadr=—e fteea(- 4 )ee-eln at 


Qo 


Qo 


1 


Aus 0, = os (01 — @o) = fe [p (G1,50) = bo| folgt daher nach (1.47) 
0 0 


1 a@ 2x 
(1.57) Spot = Soa ( agit a®D ») 
Zur Unterscheidung von der Gesamtenergiedichte (1.54) wollen wir im folgenden 
die Summe der potentiellen Energiedichte (1.56) und der kinetischen Energie- 


dichte 4? als akustische Energiedichte ¢ bezeichnen. In der linearisierten Theorie 
gilt nach (1.46) und (4.57) 


(1.58) e= “|r + (Fp +40— 2). 


200 


* Hierbei wird angenommen, daB s(r, ¢)—s9(z) und damit auch 


PQ, s, t)) —P(O1, So (t)) =P (2, t)—P(Q1, So) 


die Gr6éBenordnung O(o?) besitzt. Diese Annahme kann bei stationaren Wellen- 
feldern geniigend hoher Frequenz mit Hilfe von (1.28) bzw. (1.29) begriindet werden. 
Arch. Rational Mech. Anal., Vol. 6 17 
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Aus (1.45) und (4.58) folgt 


8s <|VOVG,+ 5 (P,+4a0—¥) (P+ 4%, —¥)| 
20 


OV, +a@— F) + QV (VO) (P+ P)—a(V®)2+VOV¥), 


1 
Sag: 

00 
also 


a ai 
(1.59) =v (@ ,+ab—P)7|_ “(VO + VOVY. 


Es sei G ein endliches regulares Gebiet mit der Randflache F. Fiir die zur Zeit ¢ 
im Inneren von G befindliche akustische Energie FE, (t) gilt nach (1.59) 


= f colt) Gee. av 
G 


(1.60) 
1 OD Gre f *aVo\ed 1 povwav. 
= [+ (@,+49 oe af—W ) leva V 


Im Fall, daB durch F keine Fliissigkeit transportiert wird, folgt aus (1.60) wegen 
aDjan =0 
(1.61) SEs) =—a [ —Wopav + f + verwar. 
dt 20 Qo 
G G 
Treten also keine 4uBeren Krafte und keine Reibungskrafte auf, so erweist sich 
E,(t) als konstant. Fiir a=+-0 kann 


(4.62) Bas (VD)? 


als spezifischer Energieverlust durch Reibung interpretiert werden. Entsprechend 
stellt 


(1.63) = VOY 


fiir VY=-0 die spezifische Energieproduktion der auBeren Krafte dar. Uber- 
tragen wir diese Deutung auf den Fall, daB 0@/én nicht fiir alle x auf F ver- 
schwindet, so folgt aus (1.60), daB der in das AuBere von G gerichtete akustische 
EnergiefluB durch 


Lurie a@ 
(1.64) Sp (t) fe (Y— @,—a®) 2° aF 
F 
dargestellt wird. Der Vektor 
(1.65) S= ~ (Y— @,—a)V@ 
<0 


kann daher — entsprechend dem Poyntingschen Vektor bei elektromagnetischen 
Schwingungsvorgangen -- als akustischer Energiestrémungsvektor interpretiert 
werden. Er gentigt nach (1.59) der Gleichung 


(1.66) Go = shy cee ae VP — * (Vo), 
Qo 
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§ 2. Stationdre akustische Schwingungsvorgadnge 


Im folgenden setzen wir voraus, daB die 4uBeren Krafte und die an den 
Randflachen fester Kérper vorgegebenen Werte fiir die Normalgeschwindigkeit 


vn zeitharmonisch mit der Frequenz@ sind. Fiir das Kraftdichtepotential Y 
kann also 


0 : 
(2.1) = + 25 Wir, t) =f, (x) coswt+ f(z) sinwt= Ref (x) en ee 
mit f=/,+7/, gesetzt werden. Hierdurch geht (1.45) in 
1 1 G2@ oD i 
(2.2) 0 Bee VO) = + + Reffe toy 
uber. 


In der vorliegenden Arbeit werden Lésungen von (2.2) untersucht, die die 
Gestalt 


(2.3) P(r, t) = U(r) coswt + U(x) sinwt= Re {U(z) eae 
mit U=U,+7U, besitzen. Die Lésungen der Gestalt (2.3) charakterisieren die 


eingeschwungenen “‘stationaren“* akustischen Schwingungsvorgange. 
Nach (2.2) und (2.3) gilt fiir alle ¢ 


1 1 02@D a o@ —iw 
eoV(5 Vo) ee ee OY 


Ce (e 


= Re{|eo7(—- VU) | o* tao / age 


2 


= cost: Relay 7(} VU) = (wm +74) U—i+ 


c2 


+ sin t-Im|957( } VU) +2 (w+ ia) Ses = 


Fir ¢=0 bzw. t=2a/2 folgt hieraus, daB Real- und Imaginarteil von 


eo7(5-VU) +S lo +ia)U—} 


verschwinden. Somit gilt 


(2.4) 0 hire VU) Uf} 
mit 
(2.5) = (w+ ia). 


Man erhalt also sémtliche stationaren Lésungen von (2.2), indem man alle 
Lésungen U von (2.4) bestimmt und 


P(x, t) = Re{U(z) e*°"} 
setzt. Daw, a und c nicht negativ sind, gilt OS arc(x?)< $a. Wir kénnen daher 
x so wahlen, da8 
(2.6) OSarex(rt)<in 
ist. 


7s 
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Beachten wir, daB die auBeren Krafte wegen (2.1) das Kraftdichtepotential 


(2.7) P(r, ) =— £ Refi f(z) e749 


besitzen, so folgt aus (1.46) und (1.47), daB die fiir das Schallfeld charakteristischen 
Gr6éBen bv, p durch 


(2.8) ee > Re{VU Cas 
und os 
(2.9) p —bo=—Rel[(a—ia) U +4 = fleio 


dargestellt werden. 

Bei dem Durchgang durch die Grenzflache F zweier Medien M, und M, sind 
die Normalkomponente der Geschwindigkeit und der Druck stetig. Hieraus 
ergeben sich nach (2.8) und (2.9) fiir rc die Bedingungen 


1 oU, 1 8U, 


Ae) est. ON ~— Qyg ON 
und aj 
(2.11) (4, — 1) U, +i f= (a, p— to) Uy +i fF. 


Hierbei bedeutet 0); die mittlere Dichte, c; die Schallgeschwindigkeit und a, die 
Reibungskonstante des Mediums M,; (¢=1, 2), wahrend U,, 0U,/én,... in der 
gleichen Weise wie ®,, 0®,/0n,... in (4.52) erklart werden kénnen. 

Fir die Punkte der Randflache eines festen Kérpers erhalten wir nach (1.53), 
falls wir 


(2.12) g(t, t) = Re{g(z) e*°4} 
setzen, die Randbedingung 


(2.13) 5 Ule) = elt). 


Wir wollen jetzt die Energieverteilung bei stationairen akustischen Schwingungs- 
vorgangen diskutieren. Wir beschranken uns auf die Untersuchung von Raum- 
teilen, in denen keine auBeren Krafte wirken. Bei stationdren Vorgangen inter- 
essieren vor allem die zeitlichen Mittel @(r) und S (sr t) der akustischen Energie- 
dichte «(z,?) und des akustischen Energiestromvektors ©(z, ¢): 


(2.14) E(t) =>— ihe (x, t) dt 
und i 
(2.15) Sy-2 fet i) dt. 
Nach (1.58) und (2.3) ist 

(=a |VOP+ 3 (@+aoy| 


oe cals, {[Re(PU eio)}2 + [Re ((a —ia) tee) 
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Aus 
[Re(A ett 2— 1/4 en tob A et)? = 1 (AreBiot 4 42 enti +4]A|? 
folgt 


(2.16) E(t) = 


4 ; 24 q2 
aq |7UP+ S$" up]. 


Im Fall a=0 spezialisiert sich (2.16) wegen x=ol]c zu 

(2.17) é(r) = i [|WU|2 + 2|U]?]. 
0 

Entsprechend erhalt man wegen 


S(t.) =——- (Gp +49)7o 


== tg [Re ((a — iw) Ue-i)] . [Re(VU ei] 


und 
Re(Ac. ()\o Reh em) == 1A Be Ve" AB hed B+ AB 
S(xr) =— ae {(a—iw) UVU+ (ati) UVU} 
0 
oe ae {a(UVU + UVU) —iw(UVU — UVv)} 
0 
= 3a, @ Im(UVU) —aRe(UVU)}, 
also 
(2.18) S (x) nate Im{(o — ia) UVU}. 
Im Fall a=0 spezialisiert sich (2.18) zu 
(2.19) S(r) = @ Im(UVU). 
20 


Das zeitliche Mittel Se des Energieflusses S;(¢) durch die geschlossene Flache F 
berechnet sich nach (1.64) und (2.18) zu 


ey hl ay 1 = 0U 
(2.20) Sz= > im { 1a) Lo Parl. 


Damit haben wir die fiir Energiebetrachtungen bei stationaren akustischen Vor- 
gangen wesentlichen GréBen @, G und S; durch U ausgedriickt. 


§ 3. Grundlésungen der Schwingungsgleichung 


Im folgenden setzen wir voraus, daB fiir die von uns betrachteten Medien 
die Dichte 9) im Ruhezustand und die Schallgeschwindigkeit c fiir |r| >R 
konstant ist und da8 die Kraftdichteverteilung / fiir |x| > verschwindet. 
Fir |r| >R geht dann die Gleichung 


(3-1) oo V(LVU) +2U =f 
Qo 

in die Helmholtzsche Schwingungsgleichung 

(3.2) AU+22U=0 
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mit konstantem x tiber. In diesem Abschnitt sollen einige spezielle Lésungen 
von (3.2) untersucht werden, aus denen wir spater die Losungen der in der Ein- 
leitung formulierten Probleme aufbauen werden. 

Wir fragen zunachst nach Lésungen von (3.2), die nur von dem Abstand 
|x —z’| =r von einem festen Punkt ry’ abhangen. Fir U(r)=y(7) erhalten wir 
die Gleichung 


” 2 , ppo ae 
(3.3) y +a tay = 0 
mit der allgemeinen Lésung 
i oc eer e tH 
(3.4) y(n) =A 4 BS 


Der Lésung e’”"/v von (3.2) entspricht nach (2.3) die Lésung 


(3.5) by 4 hem Re{ é aI — 1 6" cos (%,7 — ot) 


Y 


der homogenen Wellengleichung 


ath GON eae Oe 
(3.6) A® oh Oe fet Ore 


Wir hatten x so gewahlt, daB x, = Re(x) >0 und x,=Im(zx)=0 ist. Die Funktion 
(3.5) beschreibt daher einen vom Punkt x’ ausgehenden Wellenvorgang mit der 


Amplitude a Cat 
Die Flachen konstanter Phase x,7—qwt sind die Kugeln r= = t+ const mit 


“a 
dem Mittelpunkt x’. Ihre Radien vergréBern sich fiir wachsendes ¢ mit der 
konstanten Geschwindigkeit 


(3.7) qlee 


Der durch (3.5) beschriebene Schwingungsvorgang stellt ein stationaires Wellen- 
feld dar, das durch eine im Punkt x’ wirkende periodische Stérung mit der 
Frequenz  erzeugt wird. Fiir a=0 folgt aus (3.7) cy=c, so daB die in (1.49) 
eingefithrte GréBe c als Schallgeschwindigkeit des betreffenden Mediums inter- 
pretiert werden kann. Fiir a =£0 gilt 


(o+7a)§= + 75 Alo? + a)? + w]§ + 7[(w? + a2)4 w}*}, 
also wegen (2.5) 


(3.8) oe [oo (w? + a?) + w?]t = + £@ 


cV2 cV2 [wo (wo? a2)? — wh 


und somit 


(3.9) =o" V2 (1 +(2) 1 


vt 
2 


Es gilt also ¢y<c. Der Quotient c,/c hangt nur von a/w ab und strebt fiir a/a +0 
gegen 1, wie sich durch Reihenentwicklung ergibt. Bei auftretender Dampfung 


erweist sich somit c als asymptotischer Wert der Phasengeschwindigkeit c, fiir 
hohe Frequenzen. 
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Zur Untersuchung des mit dem Wellenfeld (3.5) verbundenen Energie- 
transportes berechnen wir den mittleren akustischen EnergiefluB durch die 


Kugel X, um yx’ mit dem Radius R. Mit Ue) =— e’*” erhalten wir nach (2.20) 


und (2.5) 
et = Im fio ia [opal 
; lt—v’|=R 
DXRE —2xRT { oF; sy —=)} 
% é m4(@ ia) (ix ; 
also 
(3.10) Srp = om en mR [o Hy + Ax%y + =) 
0 


Der mittlere akustische EnergiefluB durch 2’p ist also fiir alle R positiv, so daB 
der Energietransport zum Unendlichen hin gerichtet ist. Im Fall a=0 ist 


~ 2 
(3.11) page eels On ee 
0 Q0¢ 


Sys, ist also von R unabhangig; der Punkt r’ ist eine akustische Energiequelle 
27 w? 


der mittleren Starke 
Q0 © 


Entsprechend wird durch - e'**” baw. 


et (xr+ot) 


(3.42) D(z, t) = Re { be : é*"cos(%,7 +a 1) 


fe 


ein Wellenvorgang beschrieben, bei dem sich die Kugeln konstanter Phase 
x,7-+ot fiir wachsendes ¢ mit der radialen Geschwindigkeit —° auf den Punkt x’ 
x. 


1 
zusammenziehen. Der mittlere akustische Energieflu8 durch die Kugel 2 betragt 


~ 


INGA ° ‘ a 
(3.13) mee rr | W %,— ax, - aE 


Es gibt also ein Ry, so dab Sie fiir alle R> Ry negativ ist. Das bedeutet, daB 
bei dem Schwingungsvorgang (3.12) fortwahrend akustische Energie aus dem 
Unendlichen eingestrahlt wird. Ein solcher Vorgang laBt sich jedoch physikalisch 
durch Anordnungen, die im Endlichen getroffen werden, nicht realisieren. 

Hiermit werden wir auf eine charakteristische, von SOMMERFELD [5] erkannte 
Schwierigkeit gefiihrt, die im wesentlichen darauf beruht, da8B physikalisch nur 
diejenigen akustischen Schwingungen sinnvoll sind, die sich im Unendlichen wie 
auslaufende (divergente) Wellenvorgange verhalten, wahrend die Schwingungs- 
gleichung (3.2) auch die im mathematischen Sinn gleichwertigen, physikalisch 
jedoch nicht zu realisierenden vom Unendlichen her einstrahlenden (konver- 
genten) Wellenvorgainge umfaBt. Es ergibt sich daher die Notwendigkeit, die 
physikalisch realisierbaren Lésungen von (3.2) durch eine zusatzliche, sich auf 
das Verhalten im Unendlichen beziehende Bedingung zu kennzeichnen und die- 
jenigen Lésungen auszuschlieBen, die konvergenten Wellen entsprechen oder 
durch Uberlagerung von konvergenten und divergenten Wellen entstehen. 
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Nach SOMMERFELD kann man sich hierzu folgender Bedingung bedienen: 


Sommerfeldsche Ausstrahlungsbedingung. Die Lésung U(x) von (3.2) erfiille 
fiir |] =7—>co und |x)|=1 die asymptotischen Relationen 


(3.14) U(r tp) =0(7} 
und 
(3.15) 2 U(r te) — ix Ure) = o(4) 


gleichmaBig fiir alle Richtungen 7. 

Von den rotationssymmetrischen Lésungen (3.4) der Schwingungsgleichung 
erfiillen fiir r’=0 genau die Funktionen der Form A on die Ausstrahlungs- 
bedingung, also gerade diejenigen Loésungen, die Ausstrahlungsvorgange be- 
schreiben. 

Wir zeigen jetzt, daB die Funktion 

ei*le—-e'| 
(3.16) U(t) = eed 
auch fiir r’+-0 die Ausstrahlungsbedingung erfiillt. Bei festem ry’ strebt mit |r| 
auch |x —r’| gegen ov, so daB (3.14) gilt. Aus 


ix |e—y'| aa ts ale LH“ 1 
Aly) Vie — ei%lt zi. ) , 
ae ec tage a ee eae (e—F) 


folgt fiir 700 


=e ee ge. 
(3.18) VU( %) = ix Urn) ee, +0(5). 
Wegen 
AQ penny) olaghery ts 
ee raat aw F) 
und (3.14) gilt daher fiir 7 +o 
(3.20) V U(r tq) = 1% U(r t) vo +O (5). 


Wegen 0/¢7=r)V erhalt man hieraus (3.15) durch skalare Multiplikation mit toe 
Fir |z'| <<C gelten die angegebenen Abschatzungen gleichmaBig in Y und x’. 
Weitere Losungen der Schwingungsgleichung (3.2), die der Ausstrahlungs- 


bedingung gentigen, kann man aus den Lésungen (3.16) durch Superposition 
gewinnen. Es entstehen Lésungen der Form 


(3.21) Sy aad 
pa * [ey 


Unter ihnen befinden sich die Lésungen 


20 Ue (gy a A fem ere meal etxlens'l 
mai oi Lins le ae 


Hierbei sei n ein Einheitsvektor. Durch U, (x) bzw. 


(3.23) O: (e545 Re{U, (x) - eae. 
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wird ein Wellenfeld beschrieben, das durch zwei periodische Stérungen gleicher 
Starke, aber mit um z gegeneinander verschobenen Phasen erzeugt wird. Fiir 
t—0 erhalt man die Funktion 

(3.24) Curae eee ee 


ny’ |r—r’ ibs 


die fiir aon) wegen der Vertauschbarkeit der partiellen Ableitungen ebenfalls 
eine Lésung von (3.2) ist. U)(z) kann als Doppelquelle im Punkt r’ mit der 
Achse n interpretiert werden. Wir zeigen, daB auch U,(r) die Ausstrahlungs- 
bedingung erfillt. Hierzu wahlen wir ein a ge dessen x-Achse 
die gleiche Richtung wie n hat. Dann gilt mit r={vx, y, z} und r’={4’, y’, 2’} 


0 etxle—v'| oe ee ae 14 1 
Gets net’ a elt zi. ed ee 7 Gxt 
0 (x) Ox lr—r'| lx—z'/? \z—r'|3 ( ) 


; pele rally ea. il 
Sn 4 pe iele fay O(4)- 


Ferner folgt aus (3.17) durch Differentiation fiir 7—>co 
i) tx|y—v’| 
VUy(t) = 55, (V4 
0 (x) Ox’ t pare 


= (ix)re*le-# Tae res ola :): 


(3.25) 


(3.26) 


Aus (3.25), (3.26) und (3.19) folgt 
(3.27) V Up? t) = 1% Ulr 2) to +0 (5). 


Multipliziert man (3.27) skalar mit yy, so erhalt man (3.15). Die Giiltigkeit von 
(3.14) folgt unmittelbar aus (3.25). Fiir |x! | < C gelten die Abschatzungen gleich- 
maBig in yy und x’. _Insgesamt haben wir bewiesen: 


Lemma 1. Fiir |x'| < C erfiillen die Funktionen 


Ug)= os und Ge) = 5, Ul) 
die Bedingungen (3.14) und (3.15) gleichmafig in Yq und x’. 

Mit Hilfe der Lésungen (3.16) und (3.24) kénnen weitere Lésungen der 
Schwingungsgleichung (3.2) durch Integrationsprozesse gewonnen werden. Ist 
etwa F eine stetig differenzierbare regulare Flache und sind 4 und yw auf F defi- 
nierte stetige Funktionen, so sind 


et*le— ely 
28 y= f 40) 
oes z= mathe 
und 
ei*le—-v'| 
5 eee f Te, 
3,29) y= fae ep teeee: [leant 


fiir y¢F Losungen von (3.2). V(r) und W(x) beschreiben Wellenfelder, die durch 
eine kontinuierliche Verteilung von einfachen und doppelten Quellen auf I’ 
erzeugt werden. Nach Lemma 1 erfiillen V(x) und W(x) die Ausstrahlungs- 
bedingung. 
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Fiir die physikalische Interpretation der Ausstrahlungsbedingung ist es von 
Bedeutung, daB in gewissem Sinn auch die Umkehrung der letzten Aussage gilt. 
Hierzu beweisen wir 

Lemma 2. Es sei F eine stetig differenzierbare, reguldre, geschlossene Flache 
mit dem Aufenraum G, und U(x) eine die Ausstrahlungsbedingung erfiillende, in 
G,+F stetig differenzierbare Lésung der Gleichung AU+x?U=0. Dann gilt fiir 
reG, 


mit 


Zum Beweis betrachten wir einen Punkt x aus G,. Es sei o, eine zu F punkt- 
fremde Kugel um y mit dem Radius g und 2, eine o, und F umfassende Kugel 
um den Nullpunkt mit dem Radius 7. Ist G das von o,, 2, und F eingeschlossene 
Gebiet, so gilt 


fey aye CAO UAB ae 


n 
2,—F—o9 


Aus der Ausstrahlungsbedingung, die von U nach Voraussetzung und von Q 
nach Lemma 1 erfiillt wird, folgt 


lim (2 a u ££ \dk =0, 


eS or or z 
also tg 
ieee @ ; aQ @U ss 
[(03, UGE) ah =lim [ (048 —0 ak =— 200) 
%9 


Damit ist Lemma 2 bewiesen. 


Lemma 2 besagt, daB sich jede Funktion U, die im AuBeren von F der Schwin- 
gungsgleichung (3.2) und der Ausstrahlungsbedingung geniigt, als Summe zweier 
Flachenpotentiale (3.28) und (3.29) darstellen laBt. 


Wir wollen noch einige einfache Folgerungen aus Lemma 2 anfithren und 
beweisen hierzu 


Lemma 3. Die Funktion U(x) geniige auBerhalb der zweimal stetig differenzier- 
baren, reguldren, geschlossenen Fliche F der Gleichung AU+x%®U=0 und der Aus- 


strahlungsbedingung. Dann gelten fiir festes a und r—>oo die asymptotischen Rela- 
tionen 


(3.30) U(r t)— a) =O (+) 
039 fOr 0) x00 a) =0(4), 
(3.32) VU (1%) —ix%) U(r t,) = 0 & 


gleichmafig in Xp. 
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Aus (3.30) und (3.31) folgt, daB die Ausstrahlungsbedingung unabhangig von 
der Wahl] des Nullpunktes ist. (3.32) besagt, daB sich die Radialkomponente von 
VU fiir roo wie O(1/r) verhalt, wahrend die tibrigen Komponenten von VU 
mit der Ordnung 0 (1/7) verschwinden. Im Sinn der physikalischen Interpretation 
von VU bzw. V@ folgt hieraus nach (1.46), daB die Fliissigkeitsteilchen fiir groBes 
y nahezu in radialer Richtung zum Stérzentrum schwingen. Zum Beweis von 
Lemma 3 beachte man, daB die Flachenpotentiale (3.28) und (3.29) nach Lemma 1 
sowie (3.20) und (3.27) die in Lemma 3 formulierten asymptotischen Relationen 
erfiillen. Aus Lemma 2 folgt daher unmittelbar die Behauptung von Lemma 3. 

Als weitere Folgerung aus Lemma 2 sei noch bemerkt, daB fiir komplexes x 
mit Im(x)>0 die Ausstrahlungsbedingung mit der Bedingung U(r x9) =o (4/7) 
fiir 7—> oo aquivalent ist. 

Neben den Flachenpotentialen (3.28) und (3.29) kénnen auch réumliche Ver- 
teilungen von Punktquellen untersucht werden. Es sei G ein endliches regulares 
Gebiet mit der Randflache F. Dann ist 


rn erhle—e'| 
(3.33) T(e) = f ee’) "avy 
[g—z'| 
G 
fiirr ¢G-+F eine Lésung von (3.2), die nach Lemma 1 die Ausstrahlungsbedingung 
erfiillt. Es gilt, wie analog zu dem entsprechenden Satz der Potentialtheorie* 
gezeigt werden kann, 


Lemma 4. Es sei G ein endliches regulires Gebiet mit der Randflache F und 
t(x) eine in G+F hélderstetige Funktion. Dann ist die durch (3.3) defimerte 
Funktion T(x) diberall stetig differenzierbar, und es gilt fiir CG 


AT++2T=—4nt. 


§ 4. Die Eindeutigkeitssatze 

Die in den ersten drei Abschnitten dieser Arbeit durchgefiihrten Unter- 
suchungen erméglichen es, die Grundaufgaben der Theorie stationarer akustischer 
Wellenfelder mathematisch zu formulieren. So fiihrt die Berechnung eines durch 
eine zeitharmonische konservative Kraftdichteverteilung in einem unbegrenzten 
Medium erzeugten stationaren Wellenfelds auf das in der Einleitung formulierte 
Problem (A). Entsprechend ergeben sich die Probleme (B) und (C), falls in das 
Medium feste Kérper mit der Randflache F eingebettet sind, die zeitharmonische 
elastische Schwingungen vollfiihren, oder falls zwei verschiedene Medien mit der 
Grenzflache /’ vorliegen. 

Wir wollen in diesem Abschnitt zeigen, daB die Probleme (A), (B) und (C) 
héchstens eine Lésung besitzen. Zur Vereinheitlichung der Beweise gehen wir 
von folgender Anordnung aus: 

Die zu untersuchenden stationaren Wellenfelder denken wir uns durch eine 
raumliche Verteilung zeitharmonischer 4uBerer Krafte mit dem Kraftdichte- 
potential 
(4.1) (r,t) = — < Refi f(z)e*"} 


*x Siehe etwa [6], S.150ff. 
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sowie durch endlich viele feste Kérper Ky, ..., K,, mit den stetig differenzierbaren, 
regularen Randflachen /;, ..., j, erzeugt, die zeitharmonische elastische Schwin- 
gungen mit der Frequenz  vollfithren. Mathematisch driickt sich der EinfluB 
der Koérper K, durch die Vorgabe einer stetigen Normalgeschwindigkeitsverteilung 
der Form 


(4.2) oe Re {g(t) ef 


auf den Flachen F aus. Zur Abkiirzung setzen wir K+---+ £,=F. Das AuBen- 
gebiet G von F werde durch eine zu F punktfremde stetig differenzierbare, regulare, 
geschlossene Flache Ff’ in das AuBere G, und das Innere G; zerlegt. Den in G, 
liegenden Teil von F bezeichnen wir mit F, und den in G; liegenden Teil mit &. 
G, werde durch ein Medium M, mit der Dichte @,;(z) im Ruhezustand, der Schall- 
geschwindigkeit c; (x) und der Dampfungskonstanten a; und G, durch ein Medium 
M, mit den entsprechenden Gr6Ben @,(r), c,(z) und a, ausgefiillt. 

Wir setzen voraus, daB die Funktionen @y;, 9, und c;, c, in ihren Definitions- 
bereichen G;+ + F’ und G,+F&+F’ stetig differenzierbar bzw. stetig sind und 
daB oy, und ¢, fiir |r| >R konstant sind. Die Berechnung des durch die Kraft- 
dichteverteilung / und die elastischen Schwingungen der Kérper K; in den 
Medien M, und M, erzeugten stationdaren Wellenfelds fithrt zu folgender Aufgabe: 


(D) Es ist eine Funktion U mit folgenden Eigenschaften zu bestimmen: 
a) U erfiillt in G; und in G, die Gleichung 


Oo ies VU) +xU=f 


mit 
_J@oitt) fir xeG, | | a VeT ae 
Qo (t) = x x(r) = 
Qoa(t) fir EG, La@ria,) fix eC 


Ca(z)® 
b) U geniigt fir 7->co der Ausstrahlungsbedingung (3.14), (3.45). 
c) U ist in G;+ #+F’ und in G,+ F.+F’ stetig differenzierbar. 
d) Fir rCF’ gilt 


2 


2 “ 
(a, — i) U, +12 f= (a,— iw) Ui 


1 0U, Lita), 
Qoa On Ooi On ~ 


e) Fir r€F gilt 
0U 
on 
Die in der Einleitung formulierten Probleme (A), (B) und (C) gehen aus (D) 
durch Spezialisierung hervor: (B) erhalt man, indem man sich auf Funktionen Q 
und ¢ beschraénkt, die im gesamten AuBeren G von F unter Einschlu8 von F’ 
stetig differenzierbar bzw. stetig sind. (A) und (C) entstehen aus (B) und (D) — 
fiir nm==0 (n= Anzahl der in das Medium eingebetteten festen Korper). Die 
Findeutigkeitssatze fiir die Probleme (A), (B) und (C) folgen somit aus | 
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Satz 1. Es gibt héchstens eine Funktion U mit den unter (D) formulierten 
Exgenschaften a) bts e). 


Zum Beweis nehmen wir an, daB U, und U, Lésungen des Problems (D) sind. 
Die Funktion U=U, — U, erfiillt dann die Bedingungen b) und c) sowie 


a’) goV (2 VU) +2U =0 fiir xre€G, und réG,, 


(a, —to) U,= (a4,—i) U, 


t 


d’) {Mou Der | fir rer’, 


e’) = == 0) fir xEF. 


2, sei die Kugel um den Nullpunkt mit dem Radius 7. Wir wahlen 7 so groB, 
daB 2, die Flachen F und F’ umfaBt. Wir betrachten den Ausdruck 


(4.3) 1m {le a a) he oS art, 


der nach (2.20) bis auf den reellen Faktor $ den mittleren akustischen Energie- 
flu8 des Wellenfeldes U durch die Kugel 2, darstellt. Beachten wir, daB @, 
und » fiir |x| > R konstant sind, und setzen wir 


(4.4) Got) =e, #(e)=x* fir |x] >R, 
so folgt aus der Ausstrahlungsbedingung (3.14), (3.15) fiir yoo 


(4.5) Im {e = ss lly ie ar| = ae Tn fe*(a,+ pies o(1). 
lias t|=r 


Nach (2.6) gilt OS are x*<a/4. Hieraus und aus w>0 folgt 


(4.6) Im {x*(a, + 1o)}> 0) 
und somit fiir 7 > co 
ee ees 
(4.7) Im {(o “= eM me Pree PCW 
Si —7. 


Andererseits gilt, falls wir das von 2, & und fF’ eingeschlossene Teilgebiet von 
G, mit G, bezeichnen, 


Aus a’) folgt fiir r€G, und r€G; 


(4.9) V(— OVU) = —-|vu|?— =" |U I. 
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Das Integral iiber F, verschwindet nach e’). Fiir yc F" gilt nach d’) 
MY a 8h _, 4 2a 
(4.10) (@ — 1 @,) a U, oa (w —ta;) —U 5 


Insgesamt erhalten wir damit, falls wir noch (2.5) beriicksichtigen, 


(0 — ia) [ 7 py FM —ia,) { —|\VUP av — 
20 Qo 


t 
n 


(4.11) ler Gr gi 
4 . inane i 
—o(or+ai) f *,|UPaV + (@ sa) f Be U; 7 ak. 
ce F’ 
Aus a’), e’) und (2.5) folgt 
Uy A 
si. i Sas pgs F 
(o—ia) [ 7 OF 3 ak =(o ia) [0(- vu) V 
G; 
=(o—ia) [|-\vup—=l|up|av 
20 20 
G; 
= (o—ia) [ + |\7oPav—o(w*+ ad) [ 1_|u|zav 
20 20% 
i Gj 
bzw. durch Ubergang zum Imaginarteil 
; =a GUS 1 : 
4.12 Im !(w — ia, / Ee a. | A 
(4.12) {o Hiei |-—< rae 


Hieraus und aus (4.11) folgt 


(4.13) Im io ia) f ees aF| — an [L\pupav—a, [+ \rupar. 
iti ok or “ 20 z 20 


Da die Dampfungskonstanten a, und a; nicht negativ sind, gilt 


(4.14) Im {o Ga. f= op ar| <0, 
J @o or 
ltl=r 
also wegen (4.7) fiir r—>co 
(4.15) Im io — ia) f epee ar| =0(1) 
Qo or 
ltl=r 
Hieraus und aus (4.13) folgt 
(4.16) a, [ \Puev+a, f +\yupav=o. 
5 0 Qo 
Gq Gi 


Im Fall, daB8 a, und a; von 0 verschieden sind, gilt nach (4.16) fiir alle y aus 
G, und G; VU=0, also auch 0, Piss VU) = 0 und somit nach a’) U=0. Damit 
- . 9 

ist Satz 1 im Fall a, == 0 und a; = 0 bewiesen. 


Zum Beweis von Satz 1 fiir a,= 0 benutzen wir zwei Ergebnisse von RELLICH 
[13] und Herwz [4]: 
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Lemma 5. Fiir |x| >R gelte AU+ x2U=0 mit konstantem, positivem x. Gilt 
dann fiir roo 
f [Gilead (4) 
Igl=r 
so verschwindet U fiir |x| >R. 

Lemma 6. Es sei G ein Gebiet mit der stetig differenzierbaren, reguldren, ge- 
schlossenen Randfldche F. Ferner sei ty ein Punkt auf F und F (xy, 6) der innerhalb 
der Kugel |x —Y9| = 6 liegende Teil von F. Die Funktion U sei in G+ F (rp, 6) 
stetig differenzierbar und geniige in G der Gleichung 0,V lis VU) +x2U=0. Gilt 


dann fiir alle Punkte auf F(t), 6) U= a ==) 080 we schbanats U tdentisch in G. 


Einen elementaren Beweis des zundchst von RELLIcH [13] mit Hilfe der 
Theorie der Kugelfunktionen bewiesenen Lemmas 5 gibt MIRANKER [9]. Lemma 6 
ist in dem Eindeutigkeitssatz fiir das Cauchysche Anfangswertproblem enthalten, 
den HE1Nz [4] fiir Gleichungen der Form 


AU = F(z, iy ss 


bewies und der spater von CoRDES [2] und ARONSZAJN [Z] auf den Fall iiber- 
tragen wurde, daB an die Stelle des A-Operators ein beliebiger elliptischer Diffe- 
rentialoperator zweiter Ordnung tritt. 


Im Fall a,=0 folgt fiir yoo aus (4.15) und der Ausstrahlungsbedingung 
(4.17) fe |2 4d Hat). 
lyl=r 
Nach Lemma 5 verschwindet daher U fiir |x| >R und somit nach Lemma 6 
auch fiir rE G,. 
Insbesondere gilt auf F’ 


a) 
(4.18) Ue a 
also nach d’) auch 

ou; 
(4.19) Ce moe cg te 


Nach Lemma 6 verschwindet daher U auch in G;. Damit ist Satz 1 auch im 
Fall a,=0 bewiesen. 

Es sei noch bemerkt, daB Lemma5 nach Kato [5] auf Gleichungen der 
Gestalt 


(4.20) AU +x(r)2?U =0 


ausgedehnt werden kann, wobei x(r) fiir |z|—>co gegen eine positive Zahl x* 
konvergiert und fiir beliebiges «>0 die Abschatzung 


(4.21) (x) —x* =O (pea 


gleichmaBig fiir alle Richtungen erfiillt ist. Es ist zu vermuten, da8 sich die 
Methoden von Karo unter geeigneten, zu (4.21) analogen Voraussetzungen tiber 
0, auf die Gleichung 


(4.22) eoV(—- VU) + x(x 2U=0 
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iibertragen lassen. Es scheint daher, da in den Eindeutigkeitsuntersuchungen 
dieses Abschnitts und auch in den Existenzuntersuchungen des nachsten Ab- 
schnitts die Voraussetzung, daB x(r) und (x) fiir |r| >R konstant sind, abge- 
schwicht werden kann. Wir wollen jedoch diesen Fragenkreis in der vorliegenden 
Arbeit nicht weiter verfolgen. 


§ 5. Akustische Welienfelder in unbegrenzten Medien 
mit stetig veranderlichen Stoffeigenschaften 
In diesem Abschnitt wollen wir das erste der in der Einleitung formulierten 
Probleme behandeln. Wir betrachten ein unbegrenztes Medium mit der Dichte 
0o(x) im Ruhezustand und der Schallgeschwindigkeit c(r) und setzen voraus, 
daB 0) im ganzen Raum hdlderstetig differenzierbar und c fiir alle x hdlderstetig 
ist. Fiir |x| > seien g) und c konstant. Wir setzen 


(5.4) Go(t) = 0; C() =e Sti lit 
n 1 ; pti 
(5:2) x(x)? ae ae yk? — o as a: a) 


Wir wahlen x(r) so, daB OS arc x(x) < $a gilt. 

Die Aufgabe, das durch eine zeitharmonische Kraftdichteverteilung mit dem 
Potential (4.1) erzeugte stationdre akustische Wellenfeld zu berechnen, verlangt 
die Bestimmung einer Funktion U, die fiir alle y der Gleichung 


(5.3) eV Vu) 45? Ot 


sowie der Ausstrahlungsbedingung (3.14), (3.15) mit x=x* geniigt (Problem (A)). 
Wir setzen voraus, daB f hélderstetig ist und fiir |r| >R verschwindet. Nach 
Satz 1 besitzt Problem (A) héchstens eine stetig differenzierbare Lésung U. 

Es sei U eine hélderstetig differenzierbare Funktion mit den verlangten 
Eigenschaften. Nach (5.3) gilt 


ag 2. 4,8 1 = 
(5.4) AUS eye (x aw?) U eK as fir: - |e) de 


0 He ey eee 


Da nach Voraussetzung x und f hélderstetig sowie gy und U hdlderstetig diffe- 
renzierbar sind, ist 


(2 — 42) U QV VU Et 
0 
holderstetig. Hieraus sowie aus Lemma 4 und (5.4) folgt, daB die Funktion 


oe Ue) 
4a 
Ix’'|<R 
im gesamten Raum der Gleichung AV +x*2V =o und der Ausstrahlungsbedin- 
gung gentigt. Aus Satz 1 folgt daher V =0 bzw. 
(5.6) UG) ee fr Nea Elomi 


lv'|<R 


(x2 — 42) C= 0,V 1 VU 4 / alla av, 
Co t=" 
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mit 
(5.7) t= (x? — 22) U — o9V-VU +f. 
Da + hélderstetig ist, gilt nach Lemma 4 und (5.3) fiir |r] <R 
lott 5 U=AU +H U + (A — 0) U + eV LU 


t() + = = f {[o** Pine 


|e—z'| 
1 et * |e—2'| S hye 
Gold) * [z= ¥1 pre) ayy =f). 
Damit ist gezeigt: Ist U eine hélderstetig differenzierbare Funktion, die der 
Gleichung (5.3) und der Ausstrahlungsbedingung geniigt, so gibt es eine hélder- 
stetige Funktion t, mit deren Hilfe sich U in der Form (5.6) darstellen 1a4Bt; 
t geniigt fiir |v] <R der Integralgleichung 


see: 


— Oo (t)V 


(5.8) T(x) Tee Te) EGY) avy =F) 
mit : 

i pee ete |z—t"| 1 eee 
69) Kee) = gow] Tot el) a ee f 


Zur etierpar einer Lésung U von Gleichung (5.3) bei gegebenem oy, x und f 
diskutieren wir die Integralgleichung (5.8). Die Zuordnung 
(5.40) Be Ue Lee) K(t,v') dVy 
Talis 

kann als eine lineare Transformation im Banach-Raum C, der fiir |x| <R stetigen 
Funktionen t aufgefaBt werden*. Die Norm in Cp wird hierbei durch 
(5.11) ell = Max -0) 
erklart. Wir zeigen, daB die Transformation kK beschrankt und vollstetig ist. 
Hierzu beweisen wir 

Lemma 7. Die Funktion w(x) sei fiir |x| SR stetig. Ist ||t|]| das Maximum 
von |t(x)| am |x| SR, so gibt es eine von x und t unabhangige positive Zahl C, so 
daB fir alley mit |r| SR 


(5.12) IK r()] <Cllel 
sowie fir alle t,t, mit |t,| SR, |t.| SR und |xy,—2| <1 
(5.43) |K t(t) — K t(t2)| S$ C||t||- [41 — 22]? 
gilt. 


Beweis. Nach (5.9) gibt es eine Konstante A, so daf fir alle zx, yx’ mit 
|x| SR+2 und |z’] SR+2 


(5.14) Kerlsp op 
und ,; 
(5.15) [VK eS op 


* Zur Theorie der Banach-Raume siehe etwa [14]. 
Arch, Rational Mech. Anal., Vol. 6 48 
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Satz 2. Das in der Einleitung formulierte Problem (A) besitzt unter den zu 
Beginn dieses Abschnitts gemachten Voraussetzungen tiber Q, ¢ und f genau evne 
stetig differenzierbare Lisung U. Man erhilt U durch (5.26), wober t die eindeutig 
bestimmte Lésung der Integralgleichung (5.8), (5.9) 7st. 

Zur praktischen Berechnung von U koénnen nach Satz 2 die numerischen 
Methoden zur Auflésung Fredholmscher Integralgleichungen, etwa das Verfahren 
von CL. MULLER [10], verwendet werden. 


[7] 
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